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A TRAIN CORRIDOR SCHEDULING PROCESS 

HELD OF THE INVENTION 

This invention relates to a process for scheduling the movement of trains over a 
rail corridor having a plurality of sidings or parallel tracks with crossover switches. 

BACKGROUND OF THE INVENTION 

A rail corridor is a collection of tracks and sidings connecting two rail terminal 
areas. An example of a rail corridor 8 is shown in Figure 1, showing a single main track 
10 and three sidings 20. The western end of the rail corridor is on the left side of Figure 1 
5 and the eastern end on the right. 

Scheduling rail transportation on a rail corridor is particularly complex as 
compared to highway, water, or air transportation. Trains using a single track traveling in 
opposite directions (i.e., a meet) or trains traveling in the same direction (i.e., a pass) 
must meet in the vicinity of a siding so that one train can be sided to let the other pass. 

10 Alternatively, if there exists a double main Une with crossover switches, one train caa be 
switched to the second main line to allow the other train to pass. Also, when such meets 
or passes occur at a siding, the siding chosen must be long enough to accommodate the 
train to be sided, and the train to be sided must arrive at the siding and have sufficient 
time to pull onto the siding before the passing train arrives at the siding. 

15 The railroad must earn revenue from its transportation operations, and some of 

this revenue is generally at risk if trains cannot deliver freight on time. The destination 
time of the trains must be managed insofar as possible to prevait late penalties incurred 
by the railroad. Therefore scheduling trains across a rail corridor involves arranging 
meets and passes as required for all trains, and while also meeting the schedule for each 

20 train so that they all arrive, on time, at the end of the corridor. 
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Commercially applied scheduling processes attempted to date have been based on 
paradigms which involve simulation with branch and bound techniques to find a conflict- 
free schedule. Since a branch and bound process must sort through many binary choices 
as it proceeds toward a solution, these techniques are slow, and do not take advantage of 
5 quantitative relationships that can be adduced from the scheduling context. 

Additionally, the prior art technique search processes actually become more 
complex and take longer to arrive at a solution as the number of sidings in the rail 
corridor increases. This is due to the search algorithms that form the basis for these prior 
art techniques. More sidings requires the search algorithm to search through and consider 
10 more choices before arriving at an optimum solution. As will be shown below, the 
technique of the present invention overcomes this disadvantage. Since the present 
invention calculates a cost function where each siding represents a lower cost, having 
more sidings will make it easier for the algorithm to identify the optimal (i.e. minimal) 
cost. 

15 One prior art technique uses quantitative information such as train speed, 

destiaation, and time of departure as discrete variables in an artificial intelUgence based 
system. The artificial intelligence process involves rules that are used to search through 
the trial cases until the best case is found. In addition to the considerable time talcen by 
an artificial intelligence system to optimize a solution, it is also known that a sUght 

20 change to the initial conditions may produce a significantly different result. In any case, a 
slight change to the initial conditions will require a new and leng&y computation to find 
the optimum solution, A commercial product referred to as The Movement Planner, 
offered by GE-Harris Railway Electronics L.L.C. of Melbourne, Florida, implements such 
an artificial intelligence solution. 

25 As can be seen, tiie total set of parameters for scheduling a corridor can be large, 

and of both discrete and continuous types. Generally, a cost function based on these 
parameters can be formulated, and then some method of search is executed that will 
reduce the cost and/or find a feasible schedule for the subject trains. But, the presence of 
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discrete variables in the search space prevents or greatly complicates tibe application of 
any 'Idll-climbing" search processes based on the use of gradients 

SUMMARY OF THE INVENTION 

Cost functions that are everywhere differentiable have the advantage over prior art 
5 artificial intelligence solutions of being amenable to gradient-based minimization 
algorithms that do not have to accommodate the difficulties that arise in discrete or 
partially discrete search spaces. The present invention is a process whereby a rail corridor 
and the traia schedule along tiiat corridor can be characterized by a differentiable (i.e., 
continuous) cost function, so that a search process based on differentiation may be 

1 0 applied to scheduling train activity in the corridor. 

The present invention is an analytical process for scheduling trains across a 
corridor that is driven by a cost function to be minimized, where the cost function is a 
continuous and differentiable function of the scheduling variables. The present invention 
is an improvement over the prior art cost functions that include discrete variables and thus 

15 are not differentiable everywhere. The present invention will permit the use of search 
processes relying on gradients, and as such, will converge to solutions much more quickly 
than the prior art scheduling processes involving simulation, or searching through discrete 
options. 

The corridor scheduling process of the present invention involves three steps for 
20 identification of the optimum schedule. After an acceptable differentiable cost function is 
derived, the first step is the gradient search process wherein the gradient of the 
differentiable cost function is determined. The cost function is a sum of individual 
localizer functions. For each pair of trains in the corridor that might intersect, using the 
localizer function, the intersection point is identified as having a high value if the train 
25 trajectories do not intersect near a sidmg and lower values as the uatersection point moves 
toward any siding. The gradient process may not move all intersection points precisely to 
the center of sidings dependent upon the selected threshold value and parametric values 
of the localizer function. Instead, the gradient process varies train departure times so that 
3 
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liie set of all intersection points of trains are moved nearer to sidings. The second phase 
of the process simply moves the points precisely to the centers of sidings, selects which 
train to side, and computes exact anival and departure times for the trains at the siding to 
assure the physical integrity of the meet. Jn order to center the intarsection points at 
5 sidings and side specific trains, the speeds of die mdividual trains must be modified. This 
is accomplished during the second step of the scheduling process. 

The third step maintains the proper siding relationships between any two meeting 
trains, as determined in step two, but allows the meet time to vary in an effort to assure 
that no train exceeds an upper speed limit. This final phase is again a gradient search 
1 0 process applied to all of the meet points determined in the second step. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention can be more easily understood, and the further advantages 
and uses thereof more readily apparent, when considered in view of the description of the 
15 preferred embodiments and the following figures. Identical reference characters in the 
figures refer to identical components of the invention. 
Figure 1 illustrates of a simple rail corridor; 

Figure 2 is a string diagram illustrating the corridor scheduling problem in terms of 
intersectiag lines; 

20 Figure 3 is a flow chart for the corridor scheduling process of the present invention; 
Figure 4 illustrates the basic geometry of train trajectories; 
Figure 5 is a graph of the basic sigmoid function; 
Figure 6 illustrates the use of sigmoid sums to discriminate an interval; 

Figure 7 illustrates the construction of a localizer function from sigmoid fimctions; 
25 Figure 8 illustrates an example of a localizer fimction for two sidings; 

Figures 9 A and 9B show the modification of a localizer fimction to account for corridor 
endpoints; 

Figure 10 illustrates tiie necessary geometry to achieve a balanced localizer fimction; 
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Figures 1 1 A, 1 IB, and 1 IC illustrate a technique for approximating the economic penalty 
iRmction; 

Figure 12 shows a penalty term function for early departure of a train; 
Figure 13 is an initial infeasible string graph schedule for twelve trains; 
5 Figure 14 is a string graph for trains of Figure 13 after a gradient search of the present 

invention; 

Figure 15 shows the process whereby intersection points are moved to a siding center; 

Figure 16 shows moving the first intersection point to a siding center; 

Figure 17 illustrates the process of speed adjustments to center all meets; 
10 Figures 18A and 18B through Figures 24A and 24B illustrate certain mfeasibiUties 

created by centering meets on sidings and the resolution thereof; 

Figures 19A and 19B illustrates flie two types of siding conflict; 

Figures 20A and 20B illustrate the resolution of certam siding conflicts; 

Figures 21A and 21B illustrate the "unresolvable" siding conflict; 
15 Figures 22A throu^ 22D illustrate resolution of both types of siding conflicts; 

Figures 23 A through 23 E show the cases for downward resolvable siding conflicts; 

Figures 24A and 24B show the resolution of upward-resolvable siding conflicts; 

Figure 25 illustrates train trajectories represented as broken line segments; 

Figure 26 is an evaluation of the train trajectory vector; 
20 Figure 27 shows an adjustment of train trajectory to accommodate siding delays; 

Figure 28 shows siding details for a westbound sided train; 

Figure 29 illustrates siding details for eastbound passing trains; 

Figure 30 is a complete string graph adjusted for centered meets and train sidings; and 

Figures 31 and 32 are flow charts illustrating algorithms implemented by the present 
25 invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The traditional method of graphic depiction of a train schedule for a rail corridor 
is referred to as a string graph as shown in Figure 2. This string graph represents a time- 
5 
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distance graph of train movement in the corridor depicted ia Figure 1. The horizontal 
axis represents time, (i.e., a fixed window of time) and the vertical axis represents 
distance, with tiie point at tiie origin of the graph being the western end of the corridor, 
and the point at the top being the eastern end of the corridor. The width of the graph 
5 represents the period of interest in which the trains will be scheduled. Lines on the graph 
sloping one way represent traffic in one direction across the corridor, while lines slopping 
in the opposite direction represent oppositely-directed traffic. Only the position of the 
engine is shown. The horizontal bars across the graph, bearing reference to character 20, 
correspond to the siding locations, 

10 The invention as presented herein is described in conjunction with a single-rail 

corridor with sidings. But those skilled in the art will recognize that it can be easily 
extended to multiple track main lines with cross-over switches between the main Unes. 

The essential criterion for an acceptable schedule, as expressed in terms of the 
string graph of Figure 2, is that any two train trajectories (lines) on the graph must 

15 intersect at a siding 20. If their meets are at sidings, then in addition, a choice has to be 
made as to which train to side. 

Note that, unless all of the intersecting Unes actually intersect within the sidings 
20, the schedule is infeasible. Assuming, for the nonce, that all train speeds will be 
fixed, the departure times for the trains can be adjusted in order to move the train Unes 

20 about and attempt to place aU intersection points over the sidings 20. In another 
embodiment of the present invention, it would be possible, as well, to vary train speeds, 
which would change the slopes of the train trajectory lines, in order to place intersection 
points over the sidings 20. In yet another embodiment, both speeds and departure times 
can be varied simultaneously to find a feasible meet/pass plan for the trains. 

25 The process to be described herein treats the corridor scheduUng problem as a 

geometry, problem, rather than directly as a scheduling problem, as suggested by the prior 
art. It does so by providing a mechanism by which train trajectory lines are moved under 
control of a gradient-search process based on a differentiable cost fimction in a manner 
that moves the intersection points to or close to estabUshed sidings. 
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The search process of the present invention permits variation of speeds and 
departure times, separately or jointly, and will use an everywhere differentiable cost 
function that takes on lower values as the schedule approaches feasibility. Because the 
cost function is everywhere differentiable, an iterative, gradient-search method can be 

5 applied that assures that tlie successive schedules found by the search process in fact 
converge to a conflict-free result. 

Moreover, it is possible to include, in another embodiment of the present 
invention, the constraint that a siding must be longer than a train to be sided on it. It is 
further possible to include, in yet another embodiment, the economic costs incurred by 

10 adjusting train schedules. In other embodiments, constraints on maximum train speed 
and the early departure of trains can also be considered. 

It will be ^preciated by those skilled in the art that although Figure 2 illustrates a 
situation with three sidings and three trains traveling in each direction, the technique of 
the present invention can be easily extended to any number of trains operating in each 

15 direction and any number of sidings on the rail corridor. The concepts of the present 
invention can also be extended to a rail corridor with more than one main Une and 
crossover switches between the main Une tracks. The present invention can be applied to 
any rail corridor where one train can be switched to another track when a meet or pass 
with another train occurs. 

20 The scheduling of trains must first be feasible, but in addition, there may be 

choices as to which trains to side or the order to run trains, which helps to assure that 
economic penalties will not be uicurred or, failing that, will at least be ameliorated. 

Process 30 for obtaining both schedule feasibility and economic acceptability may 
consist of a number of steps as shown in Figure 3. First, at step 31, an initial 

25 prearrangement of the trains is done, establishing their order of entry into the corridor. At 
this point, tiie train order is based solely on due times, (represented as an input to step 31 
from block 32) with no analysis as to the corridor capacity or specific departure times. At 
step 33, an initial schedule for the trains is determined; there are several numerical 
optimization techniques that may be applied here. See for example, Numerical 
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OptimiT-ation by Jorge Nacedad and Stephen J. Wright; Springer, New York 1999; ISBN 
0-387-98793-2. 

This initial schedule is input to the gradient search process, step 34, to be 
discussed below, which minimizes schedule infeasibihty. hi anotlier embodiment the 

5 gradient search process can also minimize economic penalties incurred by the raihoad for 
the late arrival of trains and give due consideration to maximum traia speeds, early 
departure times and siding lengths. The gradient search adjusts train departure times (i.e., 
the time the train enters the corridor) and/or speeds so that meets occur near sidings. The 
process 30 loops through siding choice step 38 and the conflicts decision step 36 until all 

10 train intersections are placed at or near sidings on the rail corridor by adjusting the speed 
and/or departure time (i.e., the time the train enters the corridor) of the trains traversing 
the corridor. 

The decisions made at step 38 as to which train to side for each pair of trains 
meeting at a siding may be driven by considerations of relative economic cost due to the 
15 delays created by siding one train versus another train. This siding decision process 
represents another embodiment of the present invention and will be discussed further 
below. 

Once the siding decisions are made, some of the trajectories (those for sided 
trains) on the string graph (Figure 2) will become broken lines, (representing uifeasible 
20 meets) which may cause new schedule infeasibilities for some train trajectories. At this 
point, the gradient search can again be applied, but only to the subset of subtrajectories 
that have been driven into infeasible meets. Multiple passes through the gradient search 
step 34 and siding decision process step 38 should bring the schedule to complete 
feasibility. 

25 Figure 4 characterizes the train trajectories as lines based on the ioitial departure 

tunes (moment of entry into the corridor) and the train speed, ha Figure 4, the bottom of 
the vertical axis represents the west end of the corridor, and the positive direction along 
that axis corresponds to eastbound travel. The time window of interest for travel in the 
corridor begins at time do, and the length of the corridor is denoted by X. 
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Figure 4 focuses on characterizing one eastbound train and one westbound train, 
respectively T^- and 2}, with corresponding trajectories labeled and Zy... j/, jy denote the 
speeds and dt, dj denote llie departure times of the trains Ti and 3}-. respectively. The 
departure time of a train is the time at which it enters the corridor: for an eastbound train, 
that corresponds to a point situated on the horizontal axis of Figure 4, (i.e., t = 0) and for 
a westbound train, that corresponds to a point located on the horizontal line y^L. 

Then for traux trajectory U (eastbound), we may express the relationship between 
coordinates for any point on tiie line in the form 

t-d, ' 

or 

y = s,t-Sidi (3-1). 
For train 2} (westbound), the form of trajectory Xy can be likewise expressed as 

or 

y = -Sjt + Sjdj+L (3-2). 

We can write equations of identical form for both eastbound and westbound trains by 

writing 

y = Sit-s,d, + 0,L (3-3), 

where. the speed of westbound trains by convention will be the negative of the train's 
actual speed, and 



Jo if train2;.ise 
' ~|l if train r, is w 



^^,v...«..^.iseastbound 

' '"' ' 's westbound 



This form of a Imear equation (3-3) is not the usual form dkectly in terms of slope 
and uitercept, but in this analysis train speeds and departure times will be varied and the 
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form of Equation 3-3 has tJae advantage of expressing the train trajectories explicitly in 
terms of speeds and departure times. 

The objective of the present invention is to determine the coordinates of 
intersection points yij) for pairs of train trajectories, and move these intersection 
5 points to sidings. For trains 2} and Tj, the solution for the trajectory intersection point is 
itij,yij),whsre 

and 

5, - Sj 

(ty, yy) is derived by equating equation (3-1) and (3-2) (after making the notation change 

suggested by equation (3-3)). 

This characterization of the intersection point appUes to intersections of Uke- 
15 directed or oppositely-directed trains, so that the analysis to be developed concerning 

intersection points will adjust train trajectories involving both meets and passes. 

To this point the train scheduling problem has been abstracted to a context of 

moving uitersecting lines about until all intersection points are within certain ranges (the 

siding bars 20 in Figure 2), When all intersection points that are within the rectangle 
20 (representing the corridor 8) are also within the sidmgs 20, we have obtained a feasible 

schedule. 

It is an objective to obtain a feasible schedule using a search process that 
minimizes a cost function, and in the preferred embodiment, the preferred cost function 
will have a high value if any intersection point is outside a siding bar, and a low value if 
25 and only if all intersection points are within the siding bars. Intersection points entirely 
outside the graph are not considered; the corridor and scheduling period are considered to 
be co-extensive with the graph. 
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Let a function of a single yjj with this cost fimction property be a localizer 
fimction, and construct such a localizer function using the sigmoid fimction as a basis. 
The preferred function will depend on the basic sigmoid function, which has the equation 

<^{x\oc,0)=-^ (4-1), 
1 + e 

5 

and has a graph of the form shown in Figure 5. 

The parameter /? of the sigmoid function determines a horizontal asymptote for 
the curve, and the parameter a determines how sharply the fimction rises as it crosses the 
10 y-axis. As a approaches infinity, the sigmoid curve approaches a step fimction. In the 
preferred embodiment /?= 1.0 and 0.5. 

Because the sigmoid function can transition sharply from a low to a high value, it 
is a good continuous approximation of discrete processes. Sums of sigmoids can also be 
used to determine whether or not a variable has a value in an interval. Specifically, for 
15 ■ the interval /a, 6/. define the fimction 

D{x\a,b) = aix - a;a,p)-a{x~b',a,p) (4-2). 

Based on the graph of the sigmoid as depicted in Figure 5, the graph ofD(x; a, b) 
takes the form shown in Figure 6, which shows the function D(x; a, b) (reference 
20 character 60) derived as a sum of two sigmoid fimctions 62 and 64. 

Since each of the sigmoids 62 and 64 could be made to approximate a step 
function as closely as desired, the function D(x; a, b) can be defined to very sharply 
discriminate when ;c is in the interval [a, b], and can be made to approach a pulse of 
width 6 - fl as closely as desired. 
25 Also, since the fimction D(x; a, b) (reference character 60) approaches zero as x 

becomes more distant firom the interval [a, bj, it is possible to sum such iaterval 
discriminators (for non-overlapping intervals) and thereby obtain a function which takes 
a high value when x is in any of the intervals of interest, but is low otherwise. This is 
shown in Figure 7 for the two intervals [aj, bj] and [a2, bz], and it is obvious to those 
3 0 skilled in the art that the construction is generalizable to any finite number of intervals. 

11 
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The localizer function 70 shown in Figure 7 (generated by summing sigmoid 
functions 72, 74, 76, and 78) can be extended to any fibite number of intervals, so such a 
localizer can be constructed for any corridor of the type in Figure 1 (one main track, one 
or more sidings). The sidings are represented along the x-axis between points ai and bi. 
5 The localizer function 70 has the form 

L'{x;a,p,a,,b,,a2,b2)==fi~i,tr{x-ai',a,i3)+Y,(r{x-bi',a,fi) 

1=1 M 

The cost function for the scheduling problem of Figure 2 will be derived below 

10 usiag the localizer function concept, and assuming ns sidings. Iti the preferred 
embodiment, the cost function will be low if and only if the y-coordinate yg for an 
intersection of train trajectories lies within the range of a siding, but the localizer function 
70 of Figure 7 in fact displays the opposite effect. Thus we will first define the locahzer 
function 

15 

LXx;a,fi,a,,b, «,,,&„,)- - gcrCx - ar,a,/3)+ f^cx{x - br,a,fi) (4-3), 

which has the desired property of taking a low value if and only if x is in one of the 
intervals 

20 [«p6.1-.k.^«J' 

and a high value otherwise. That is. Equation 4-3 defines a localizer function that is the 
inverse of the localizer function 70. See the localizer function 80 in Figure 8. 

The localizer function as defined above in Equation 4-3 (and taking the form of 
25 the inverse of the localizer function 70 in Figure 7) will now be used to define a cost 
function which takes lower values as the intersection points of train trajectories are 
moved toward sidings. Two versions of the cost function are described separately below. 

A Simplified Feasible-Schedule Cost Function 

30 
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"Now letting «r be the set of all trains to be run in the corridor, and lettmg Li 
represent the train trajectory line for train Ti (as in Figure 2). Define a set I of all 
possible y-coordinates of the intersection points between the train trajectories by 

1= k 'M=Li ^''j ^ {1 (^-^>- 

Note that, with reference to Figure 2, this set includes all possible intersection 
poiQts between train trajectories, even though some of those points may not be within the 
corridor 8 and/or time window of interest. It is necessary to consider such out-of-corridor 
intersection points because the search process will move the train trajectories, and may 
bring into the corridor 8 an intersection point that initially was outside the corridor 8. 

To create a cost function that takes on a low value if and only if all intersection 
points lie within one of the sidings 20, we sum localizer function values derived from 
Equation 4-3. Specifically define the vector that represents all intersection poiats in the 
vicinity of 

y = {y,ly,el ) (5-2), 
and define the cost fiinction C'(y) by 

C'0^)= Y,L{y,;a,p,a,A.-,a„^A.) (5-3). 

The cost function is a multidimensional function of the vector y, where each value 
of the vector yields a different sum based on locaHzer function values. Each localizer 
ftmction value comprising the sum indicates whether an intersection point is in a feasible 
range (the siding bars 20 of Figure 2) or not. See the cost function 80 of Figure 8, where 
the X-axis represents distance along the corridor. If all of the uitersection points relative 
to a specific siding are feasible, C'(j5) should take a low value in the vicinity of points 
represented by that sidmg; otherwise, it takes a value near tiie value of p. When many 
intersection points are involved, y(?may have to be chosen so that the near-zero sums of a 
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large number of feasible intersection points do not result in a value in the range of y5, 
which would mask the feasibility that is to be discriminated by the function. 

C'(y) is a differentiable fimction of the vector y (the intersection points) and 
therefore in each of the variables that determine the various intersection poiats, i.e., 

5 departure times and/or speeds of the trains. Therefore tlie cost function can be used with 
gradient search technique or other search techniques based on partial derivatives, to 
minimize tiie cost function value at sidings. One such technique will be discussed below. 
Since each intersection point occurring as a component of j; is a function of the train 
departure times and the speeds of the correspondiag trains, we may treat the cost ftmction 

10 as one which may be optimized by adjusting either speeds or origination times of the 
trains, or both. 

Accounting for Corridor Endpoints 

1 5 The fact that the intersection points in I may not always represent intersections, of 

trajectories within the corridor 8 poses a difficulty for the cost function as defined in 
Equation 5-3, which is that any intersection point outside the corridor is a "don't care" 
point for the search process (so long as it remains outside the corridor), but the cost 
function as defined in Equation 5-3 will assign a high value to such a point. Recall tiiat 

20 tiie cost function of Equation 5-3 is based on the localizer function of Equation 4-3, 
which is illusti-ated by reference character 90 in Figure 9A. Thus as Equation 5-3 is 
currently formulated, an otherwise feasible solution might be masked by such a "don't 
care" point. 

In another embodiment of the present invention, the solution involves modifying 
25 the localizer function 90. Figure 9A depicts the localizer function 90, as defined by 
Equation 4-3, and a modified localizer function 92, which is generated by adding two 
more sigmoid functions 94 and 96 to account for the end points of the corridor 8. 
Specifically, define 

e = y-coordinate of the eastern end of the corridor 8, 
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w = y-coordioate of the western end of the corridor 8, 
then alter the definition of flie localizer function by including the sigmoid functions 74 
and 96 as follows. 
5 L{x\a,p,a^,b^,...,a„^,\^,e,w) = L\x\a,p,a^,b^,...,a„^,b„^) 

-(T{w-x;a,p)-a{x-e\a,p) (5-5). 

The use of the localizer function L of Equation 5-5 also requires rewriting of the 
cost function in Equation 5-3, as follows. 
10 C(y)= Y,Li^iP<^,P.a„b„...,a„^,K,>e,w) (5-6). 

This cost function then should take a high value so long as any train trajectory 
intersection point within the corridor is infeasible, but has a low value for all feasible 
intersection points, as well as intersection points that fall outside of the corridor. 
15 Like C'O?), Ciy) is a differentiabie function in each component of the vector 

{y). Any Gradient search techniques or the use of other information based on partial 
derivatives, can be used to minimize the value of C{y) in the regions of the sidings. 

A Balanced Feasible-Schedule Cost Function 

20 

As can be seen from the localizer functions 70, 90, or 92, the sidings (as 
represented by the x-axis values ai to bi) are shown as being of different lengths. In fact, 
rail corridors typically have sidings of different lengths. The consequence of different 
length sidings, with respect to the cost function (see Equation (5-6)) is that the cost 

25 function minima corresponding to sidings do not have the same y value. See the cost 
function 80 of Figure 8. For siding SI, the cost function y value is represented by 
reference character 82 and the y value for siding S2 is represented by reference character 
84. Note that the minimum at reference character 82 has a larger value than the minimum 
at reference character 84. Because the sidmgs are different lengths, the sigmoid sum that 

30 is creating the minimum uses a narrower portion of the sigmoid fimction for narrower 
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sidings and, therefore, tihie associated minimum does not drop down as far as for a wider 
sidiag. This effect may cause fee cost fimction gradient optimization process to favor a 
long siding with a deeper rninimum when it is located very close to a short siding wife a 
shallow minimum. In fee embodiment discussed below, fee cost function will be 
5 adjusted to achieve equal minima for all sidings. 

If the derivative of fee localizer function 80 has a zero exactly at the midpoint 
between sidings, feen the search process will have no tendency to favor one siding over 
another. We will call such a locahzer function balanced. The situation depicted m Figure 
8 does not assure that fee cost function derivative will have a zero properly situated; 

10 alfeough fee derivative may appear to be zero between the sidiags, it can be shown by 
feose skilled in fee art, through equation manipulation, feat fee zero is usually off-center. 
Figure 10 illustrates a means to achieve a close approximation to a balanced cost 
fimction. In Figure 10, fee intervals [c;, bj], [02, b2], and [03, 65] represent locations of 
sidings along fee main corridor. We would like to assure that fee derivative of fee 

15 localizer function, as defined for this conidor, will be zero at fee midpoints mi 2, and m23 
between sidings. The localizer function generating the cost function is a sum of 
sigmoids, each of which contributes substantially only within the immediate vicinity of 
the sidings for which it creates a minimum in the locahzer function. If we assume that 
fee localizer function at point m]2 does not depend significantly on the sigmoid terms 

20 other than feose used to create minima for fee two immediately surrounding sidings, then 
we may write a simplified locahzer in the form 

(5-7). 

Note here that the sigmoid functions used to generate the locahzer function are only those 
25 sigmoid fimctions representing sidings to fee left and right of fee point of interest on fee 
localizer function. 

It can be shown by computation feat 
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provided tibiat 

-\ = a^- m,2 & Tn,^ - Cj = 63 - m^^ , 

5 as is shown in Figure 10. This requirement will force both sidings to be of the same 
length, of course, and additionally, the next siding corresponding to the interval [cj, 63] 
must then be the same lengfli as the siding corresponding to interval [a^, 62]. It follows 
by induction that all sidings along the corridor must have equal lengths for tiie locaUzer 
function for the corridor to be balanced, 
1 0 Bringing such an artifact to bear would have two effects: 

(1) the search might, at least shghtly, mislocate intersection points, since the exact 
position of the sidings would not be reflected in the model; 

(2) siding lengths would not be accurately represented relative to train lengths. 

Of these two drawbacks, the latter is in fact of no consequence, because the 
15 modification to the localizer function to account for siding lengths will not affect the 
subsequent step of the present invention (to be discussed below) wherein train lengths are 
considered relative to siding lengths. The former effect will be of minor consequence, 
since getting train intersection points nearly into the vicinity of sidings will allow minor 
adjustments to train speed to ensure intersections occur at sidings. This step of the 
20 present invention will also be discussed further below. 

In another embodiment especially favorable if there is a large discrepancy between the 
shortest and longest siding, begin with all sidings assumed equal, thereby preventing bias 
among sidings in die early part of the search, and then adjust the localizer slowly back 
toward correct siding lengths as the search process iterates. 
25 Specifically, this may be implemented in another embodiment of the present invention as 
follows. Before the search process begins, 
(1) compute the average siding length Savg as 
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'■^-^ (5-8); 

(2) redefine the position of each sidiag 5',- (corresponding to corridor interval [a,-, fe,]) as 
corresponding to the interval [a'iiO),bliO)], where 

,,^^d±:f^,^^ = ^l^±tI^ (5-9); 
2 * 2 

(3) define, for any integer n > 0, 

a;(n) = a;e-^ +a,(l-e-*') and h\{n)=b\e'^ +6,(l-e-*') (5-10), 
where A, is a positive real number. Note then that 

(3!(o) = a\ and lim a'.{n) - a, , 

and 

Z>;(0)-6; and lim &;(«) = 6,. . 

Begin the process by letting « = 0 and then as the search proceeds, increase n according to 
some scheme. 

For example, one preferred scheme would be to note when successive values of 
the cost function (during the gradient search process as discussed below) have a 
difference smaller than a predetermined threshold (see for example, the threshold value b 
referred to in conjunction with Equation 8-3 and the textual material following 
immediately thereafter), then begin to increase n (relative to the differences in siding 
lengths) and recompute the localizer function until siding lengths are within 5% of being 
accurate. This will permit the initial localizer to correspond to the balanced localizer, so 
that sidings will not tend to be favored solely by length. The initial "push" of 
iatersections toward one or another siding wiU be unbalanced. As n increases, the 
locahzer function will more accurately reflect the true corridor structure, so that 
eventually an accurate schedule is obtained. 
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Accounting for Train Lengths versus Siding Lengtlis 

The cost fimction as described above permits a search for a feasible schedule only 
insofar that trains will meet in the vicinity of sidings. No reference has been made to the 
lengths of tlie trains relative to the sidings, and if two trains have a "feasible" meet at a 
siding that will hold neither of them, then the situation is not actually feasible. There are 
other reasons that trains may not use a siding, related to grade, transportation of 
hazardous materials, etc., so the following analysis to block the use of a siding by a given 
train refers to more situations than just train length versus siding length. 

The cost function of Equation 5-6 wiU not prevent such an rnfeasibility from 
occurring, but in another embodiment, a simple modification of the localizer functions 
(Equation 5-5) on which the cost function is based will suffice to prevent such 
infeasibilities. 

In particular, the cost fimction contains a term for each possible train trajectory 
intersection point. In the previous embodiment all such terms are of exactly the same 
fortn. Now suppose that we define the locaHzer functions to be specific to each possible 
intersection point of train trajectories, as follows. la Hiis case, we generalize firom the 
context of Figure 2, and assume a total of ns sidings Si,...,S„^ along the corridor, and nr 
trains T^,...,T„ . We need the following notation: 
let 

Hi = the length of siding Si (i = l ns} (6-1), 

and let 

Mi = the length of train Tt (i =1, ....nr) (6-2). 

For any two trains 2} and 3}, define the following set of sidings from amongst all sidings 
in tiie corridor: 

Sij= te/;te {l J\^}&((M, <H,)v(Jkf^. ^H,))] (6-3). 
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Sij is the subset of sidings along the corridor on which at least one of the two trains Ti 
and Tj can be sided. Now, if the localizer function for the intersection point of the train 
trajectories of T, and 2} does not include the sigmoid terms (see Equation 5-7) 
corresponding to sidings not in 5ij, then it will remain high even though yy is within a 
5 siding, but the siding is too short for either train. In this way, accounting for siding 
versus train length actually reduces the computational complexity of the cost function. 

To specifically redefine the cost fimction in tiiis form, first redefine the localizer 
fijnctions to be specific to train pairs, i.e., 

heSf, heS„ 

10 (6-4)- 
Where the subscript "h" identifies a siding. 

Finally redefine the cost fimction as 

15 which extends the definition of feasibility so that now the value of C{y) will be low if 
and only if 

(1) all tram trajectory intersections occur on siding bars, and 

(2) at least one of the two trains in such an intersection can be sided in the corresponding 

siding. 

20 Note that this technique can be extended beyond the consideration of train length versus 
siding length: if neither of two trains Ti and Tj can be sided on siding Sk for any reason, 
then the localizer for the intersection point yy should omit the term corresponding to Sk. 
For example, we may have a case where a coal train could be sided at Sk, but would be 
unable to restart because of grade, but the interfering train, a multimodal, is absolutely not 

25 to be sided for a coal train. In this case, tiie siding may be long enough for eitiier tram, 
but would be precluded from consideration anyway. Clearly in other embodimesnts the 
defuiition of each Sy can be contracted to exclude cases such as this, thereby sharpening 
the ability of the search process to prevent vmacceptable sidings. 

20 



wo 01/49550 



PCT/USOl/00080 



Economic Costs, Early Departure, and Speed Constraints 

The cost fimctioii as described by either Equation 5-6 or 6-6 will facilitate the 
5 finding of feasible train schedules, but includes no cognizance of the other effects of 
altering individual train schedules to achieve feasibility. In another embodiment, the cost 
function is modified so that it jointly considers schedule feasibility, and the economic 
cost of late arrival. 

1 0 Economic Costs (i.e., Late Arrival) Function 

Raikoad freight service may incur various types of incentives for on-time delivery of 
freight. For the moment, consider just two types of delay penalties: 

(1) step function penalty - if a train li misses a preset delivery time tj , there is a fixed 
15 penalty cost /t// 

(2) step fimction plus linear increase - if the preset dehvety time ti is missed, there is an 
immediate penalty hi (possibly 0) which thereafter linearly increases at a rate of mt 
dollars per hour. 

Figure 1 1 A depicts a single generic form for both of these cases, since both hi and mi 
20 may be zero or positive. Thus, Figure 11 A illustrates a combined penalty function 
including both a step penalty plus a linear penalty. 

The cost function as proposed is in not a differentiable function since it lacks a 
■ defined slope at the time ti. This fact precludes, or at least complicates, use of any 
gradient search technique for ininiimzing econonoic cost unless special allowances are 
25 made at or near the tune tt. For this reason, Figures UB and IIC depict two 
approximations to the cost function, a step plus linear penalty, and a linear penalty only, 
respectively. In both figures, a line segment is grafted onto a sigmoid function in such a 
way that the resulting function remains differentiable at aU points. 

For the step plus linear penalty, a sigmoid is used to represent the cost up to a time 
30 slightly beyond tj, to which is then appended a line of slope mi. See Figure UB. 
Provided the crossover pouit from sigmoid to line segment is chosen at the point of the 
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sigmoid where the slope is exactly mi, (reference character 110) the resulting 
approximation is differentiable at all points, and therefore smoothly integratable into a 
gradient search process. If {tc, yc) represents the crossover point; then the differentiable 
version of the penalty function may be defined by 

5 

f /!,crO fort <t„ ,„ ,^ 

4(f;f„A,,m,)= ' ^ ' ^ ^ (7-1) 
^ ^ ' ' ' [m^t + y^- niit. for t > t^ 

The sigmoids used here wiU all have Pi values of 1, so that the notation for the 
parameter /?i in each sigmoid will be suppressed. This choice is made so that the 
10 asymptote of the sigmoid is determined to be A,-, in conjunction with the penalty value to 
be represented. 

The value of a,- is positive, and may be chosen to approximate the step cost as 
sharply as desired. In one embodiment the search is started with "gentle" sigmoids, then 
increase the values of the a/'s as the search progresses. This allows the early search to 
1 5 progress toward correct economic decisions rapidly, and then in the later stages of search, 
the Information concerning economic cost is sharpened to provide more accurate final 
results. 

In order to determine the crossover point 110 {{tc, in Figure IIB), it is 
necessary to solve the equation 

20 |M^;a,)| (7-2) 

for the value of t^, with tc > U. The technique for solving this equation is well known to 
those skilled in the art. It should be mentioned that the slope of <T(?-?,.;a,.) is 
everywhere positive, and takes a maximum at flie point t = U. That maximum can be 
25 driven as high as possible by selecting a large a,-, so solving Equation 7-2 is always 
possible. 
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Finally, for the purposes of expressing the gradient as will be explained below, 
note that the independent variable t in Equation 7-1 is in fact a function of the departure 
time di and ^eed j,- of train Ti, and therefore we may rewrite the equation as 



Atis,,d,;ti,h„m,) = 



{oxd.+ — < 



L 



(7-3). 



Figure IOC also iises a transition from sigmoid to line segment at point 112 on the 
sigmoid where the slope is exactly that of the Hne: the difference is that in this case the 
crossover point tc is less than U. Except for that fact, the approximating function has a 
description identical to that provided in Equations 7-1 and 7-3, 

Now, we extend the cost fixnction of Equation 5-6 or Equation 6-6 as follows. The 
extended cost function accounting for both schedule feasibility and economic cost is 
defined by 

Fiy)= vciy)^ {i-n^LA^MA,^) (7-4), 

where 

e [0,l] is a weighting factor between 0 and 1 used to adjust the relative 
importance between economic and schedule feasibility considerations. 
d = {di,d2,:.,d„^ ) is the vector of train departure times, 
s ={sj,S2,.:,s^^) is the vector of train, speeds. 
In fact, the intersection points y of train trajectories are functions of the tram departure 
times and speeds, so we may rewrite Equation 7-4 in the form 

F{s,d)= 7]C(^,d)+ (1 - T])J^Af{si,d,;t,,h,'^i) (7-5), 

and it is from tiiis latter form that the gradient may be directly computed as discussed 
below. 
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The value of the weighting factor 77 must be chosen, and the choice is of some 
importance. Note that tiie cost function as defined in Equation 7-5 will be driven upward 
both by infeasible scheduUng choices, as well as choices that make trains late, and vice 
versa. The difficulty arises when chaages of departure times or speeds cause 
5 countervailing effects in the two halves of the cost fimction of Equation 7-5. If the first 
term, representing feasibiUty, is driven up by less than the second term, representing 
timeliness, is driven down, then the search process may be emphasizing economic cost to 
such an extent that it converges on infeasible schedules. 

In one embodiment, the weighting actor 77 can be varied during the search. For 
10 example, startuig with a low value of rj would tend to try to force low economic cost at 
the expense of feasibility. This might cause the trains to swap places in the lineup, to 
improve the overall timeliness of arrivals, before the actual emphasis begins on selecting 
speeds and departure times that create a feasible schedule. Li any event, the decision as to 
how to vary 77 during the search will benefit from actual testing with examples, and final 
15 mechanism for modulating 77 will necessarily come from experience familiar to those 
skilled in the art. 

An approximate process for gaugmg the weighting factor 77 is to note that the cost 

components C(?,tf)and ^^(s;,d';;r,.,/j,.,w,) comprise different numbers of summands, 

and therefore have different magnitudes approxhnately in proportion to the number of 
20 summands involved. For example, if there are a total of twenty trains, resulting in sixty 

intersections on the string graph, thenC(5, J) comprises sixty summands and 

^A(s;,di;ti,h;,mi) comprises twenty summands. To more or less equalize the effects 
of these two contributions to the cost function, one would set the weight 77 to the value 

24 
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77 = 20/(60+20) = 0.25, thereby equalizing the contribution of each half of the cost 

function (i.e., the two terms C{s,d)md ^A{Si,di',t„hi,mi)) to the total cost. From 

this example, it can be seen that the establishment of a specific value of 7 is very specific 
to the situation under study, as is generally recognized by those proficient in the art of 
5 complex optimization. 



Early Departure Cost Function 

The late penalty assessed for economic reasons will tend to prevent train 
departures firom being arbitrarily late. However, the formulations of cost fimctions so far 
given (Equations 5-6, 6-6, 7-5) have no terms which prevent the train trajectories firom 
being arbitrarily early. A cost function to prevent early departures can be formulated in 
terms of the ubiquitous sigmoid function by defining a cost 

Ei^y%~a{d-e,;a;)] (7-6), 

where 

e,- is the earUest possible departure time for train I/, 

and 

a] is denoted with a prime to distinguish it from tlie of Equation 7-3. 
Figure 12 represents a term of this cost function for train Tc, clearly, it rapidly becomes 
high as train 7/ is pushed toward an unrealizable departure time, and rapidly drops as the 
departure time enters the realizable region. There is no actual economic cost associated 
with early departures, just a feasibility issue. Therefore the terms of Equation 7-6 
representing each train are arbitrarily given a height of 1, (i.e., sigmoid asymptotic value 
of 1) and this Equation 7-6 can likewise be combined with the cost functions for schedule 
feasibility and economic cost. Specifically, let 

G(?,rf)= 77,c(?,rf)+ (7-7), 
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where 



771+772+773=1 



(7-8). 



In one embodiment, the specific weighting of the components of cost in Equation 
7-7 can be calculated as discussed above in the example with twenty trains and sixty 
intersections points in the corridor. The schedule feasibiUty and early departure terms 
will each have sixty sxmmiands and the economic penalty term will have twenty terms. 
Using an equation similar to the one set forth above for calculating 77, we calculate 
77, =1/7, 772 =3/7 and 773 =3/7. Other weighting values can be established based on specific 
user circumstances. 

Maximum Train Speed Cost Function 

In the embodiment when the search process is permitted to vary train speeds in 
order to achieve feasibility and cost minimization, there must be a means to prevent the 

speeds from exceeding practical limits for the trains and tracks involved. In this 
embodiment we will create an additional component of the cost function that will enforce 
such speed constraints. Such a speed constraint can be implemented analogously to the 
early departure constraint of Equation 7-6. Specifically, defme a speed cost function as 



gl'^) = the maximum allowable speed for train T/. 

Like the other cost functions discussed herein, since the maximum speed cost 

function is derived from a sum of sigmoid functions, it id a differentiable function with 
respect to the intersection points of the trains on the conidor. Therefore a gradient search 
process can be used to find the minima of the cost function values. 
The total cost fimction, including feasibility of meets and passes, constraints on early 
departures and late arrivals (i.e., economic penalty), as well as constraints on maximum 
train speed then is a generalization of Equation 7-7, namely 




(7-9), 



where 



G(?,rf)= 77,c(?, J)+ %Ai^,d)+ nA^h ■n.v^) 



(7-10), 
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where 

771 + 772+773+77^=1 (7-11). 
The specific values of the weighting factors for tbs components of Equation 7-10 
can be determined by experiment. In one embodiment, using the same scheme as set 
5 forth above in conjunction with Equation (7-8), for twenty trains and sixty intersections, 
77, =0.1and 772 =773=774=0.3. 

The Gradieut Search Process 

10 The gradient V/(3c)of any function f(x) is a vector in the same space as the 

independent variable 3c which points in direction of maximum change off(x) within a 
small local area on the function's surface, thereby pointiag the way toward a local 
minimum or maximum. As such, it is much heralded in the legends and poetry of 
optimization theory. Calculation of the gradient of the various cost functions discussed 

1 5 below will permit location of the local minima identifying schedule feasibility. 

ha the current context of train scheduling, as will be appreciated by those skilled in the art, 
there are a number of possible parameters describing a train trajectory that may be varied 
to resolve conflicts within a rail corridor, i.e., to drive the cost function lower. The 
mathematics for a gradient search varying only the departure times or speeds of trains, 

20 and then varying both departure times and train speeds is discussed below. We will first 
deal only with the cost fimction associated with schedule feasibility (Equation 5-6) but 
will then extend the cost fimction to include considerations of economic costs, early 
departures and maximum train speed as discussed above, and represented by the cost 
function of Equation 7-10, 

25 

Gradient Search to Optimize Schedule Feasibility By Varying Only the Train 
Departure Times 
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First, assume that there are nr trains and permit the vector y (as represented in 
Equation (5-2)) to contain all possible intersection points. But each intersection poxatyy 
has the characterization given in Equation 3-6, which is repeated here for convenience. 

* Si - Sj 

ihsnyij is directly expressed ixi terms of departure times and speeds for all of the trains in 
the schedule. Also for convenience, recall the notation for speeds and departure times 
originally introduced above, which are repeated below. 
L = the length of the corridor, 

Si = ttie speed of train Ti (taken as a negative value for 2/ westbound), 
di = the departure time (time of entry mto corridor) of train tu and 

fo for Ti eastbound 

' [1 for T. westboimd 

Next define the vectors s = [si s„J and d = {d^,...,d„^). Express the cost 

function in tihe following terms. For notational convenience, suppress the dependence of 
the localizer and cost functions on a and p. 

C(y)=c(s,d)= I,L,(s,>SjA,dj) (8-2). 

The objective is to vary the vector d (train departure times) in order to drive the cost 
fimction lower, and one technique which can at least locate a local nunimum of the cost 
function is the gradient-directed descent, defined iteratively as follows. 

(1) Start with an initial estimate for departure time, d^ for each train «r, a stopping 
criterion, s>0, and a step size h. 

(2) For the estimate, d„ , compute Ihe gradient V(%(rf),^, of the cost function at d, , 
for varying only d , and normalize it so that it has an absolute value of 1, i.e., define 
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(8-3). 



In flie notation, the dependence of the cost function on s is suppressed, since for the 
nonce we are varying only d . 

(3) Compute the value C„=c(j„), compute d„^^=d„-hg, and then compute 

(4) If \C„ - C„+i| < e , then the search is stopped, and d^^^ is accepted as the final answer. 
Ottierwise, replace d„ with rf„^.j and return to Step (2). In the preferred embodiment, 
the search is stopped when \C„ - C„„| ^ (.001)|Co - Qj . The stopping threshold for 
such problems is very situation dependent, as is generally recognized by practitioners 
of the art of optimization. 

It remains to explicitly represent ttie gradient which is used in the 

iteration. The cost function as shown in Equation 5-6 is a fimction of the vector of 
intersection points, y , and the components of y are functions of the components of the 
vectors ? and ^ . Since at tliis point only d is variable, the gradient V(^b(rf)ofthecost 
function, is a vector of the form 

and we may obtain each component of by applying the chain rule for 

differentiation; 
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_ (8-6). 



Using Equation 8-6 and Lenuna A3 in liie Appendix A we can finally express the 
Jt-th component of the gradient as 

~-te)=^ S {S[cT(y.-^.)(^-cr(y,-6j-cr(y,-a,)(^-a(y,-«> 

c^i^--yv.JfiM^-yik))-<^iy-^ -.)(^-cr(y, (8-7). 



Gradient Search to Optimize Schedule Feasibility by Varying Only the Train 
Speeds 

Much of what was developed above can be applied here as well. The primary 
15 difference is that we now emphasize that C(y)may be regarded as a ftmction of the 
vector i , with d held constant, and we wish to vary s to seek a local minimmn of the 
cost function, and suppress the dependence on . We may therefore represaat C{y) as 
C(y)=C(5) (8-8). 

Starting with 



(8-9), 
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we need to compute the components of the form -^ic(s)) using the differentiation 
chain rule. To that end, note that 

Z^.U)Vf f Z^.(y.)V I (8-10). 

and we proceed to obtain an explicit expression for -—(yit), as follows. 



ds^^'"^ dsX s,-h ) 



W,{d-d,)^{s,e, -s,e)l-\)+{s.-s,){s.{d., ~d,)-0A 

Exploiting Equations 8-10, 8-11, and Lemma A3 of the Appendix A provides a 
final explicit form for the gradient, as shown below 

P i=l 

cT(w->.,)(^-cr(v.-j;J)-cr(y« -^)C5-cr(y, (8-12). 

The search rule using the gradient as computed in Equation 8-12 is an exact 
analogue of the search rule given in Equation 8-7 with any occurrence of the vectors 
d,d„d„,d„^, . . . 

replaced with the vectors ?,?(,, ?„,?„+i . . ., respectively. 

Gradient Search to Optimize Schedule Feasibility by Varying both Departure Times 
and Train Speeds 
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Bearing in mind that the speeds and departure times of trams can be varied 
independently, we may also exploit the expression of the cost function as a function of 
both s and d , i.e., 

C(y)=c{s,d), 

and consider the joint variation of speed and departure time to seek a cost function local 
minimum. &i this case, the gradient vector takes the form 

v(«)4-..-)=(Ac(Ml...,^4v)i^c(Ml-..^4.^)] 

Because s and d are not functionally dependent on each other, it follows that 

A^{44.±m and^(c(...1=^(c(4 (8-14). 

so the components of the gradient of Equation 8-13 are akeady determined by Equations 
8-7, 8-11, and 8-12. 

The search rule in this case is of the same form as Equation 8-7, except that we 
consider the aggregate vector 

v = {s,d) ■ (8-15), 

and replace all references to djQ,d„,d„^^va. that rule with references to v,Vo,v„,v„+i, 
respectively. 

Including Early Departure Effects in the Gradient Search 

Recall from the discussion above that a cost function causing high cost for early 
train departures, and low cost otherwise, can be posed in terms of the sigmoid fimction. 

Repeating Equation 7-6, 
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where 

et is the earliest possible departure time for train T(, 

and 

5 a'l affects the steepness of the rise of cost as early departure is 

approached. 

The value of a' may be set by experiment, but the results should not be 
particularly sensitive to its value. A good first guess, in one embodiment, for the value of 
a'l would be 0.8, although this parameter might be made smaller if there is some latitude 
10 as to the earliest departure times. 

If we wish to combine this early departure cost with schedule feasibility cost, we 
do so in a weighted sum of terms, i.e., 

D(s,d)=T]C(s,d)+il-v)E(d) r]B{0,\) (8-17), 
This concept was previously discussed above. See for example, Equation 7-7 where the 
15 aggregate cost function includes schedule feasibility, economic cost, and early departure 
effects. 

Since the gradient operation is linear on the space of functions to which it applies, 
we may write 

v(d(?,^))= 77v(c(j,J))+ (1 - ?7)v(e(^)) (8-18). 
20 We will rely on the previous gradient computations for the first term on the right 

side of Equation 8-17. See Equation 8-7 with the substitutions set forth in Equation 8-15 
and the text following. 

To deal with the second term on the right side of Equation 8-17 or 8-18, assume 
only the departure time vector d wiU be varied in a search for a schedule which is both 
25 feasible and prevents early departures. We then wish to deterraine the gradient of E{d) 
relative to llie vector d , which is of flie form 
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and (see Equations 8-6 and A-2) 

=-(i-,Xi-c'fe-^.;«:>-fc-«.;«;)) (8-20)- 

We can now construct the gradient using Equations 8-7, 8-18, and 8- 

20. Departure times are independent of train speeds, so the cost component E{d) does 
not depend on tiie speeds s . Thus the final form of the gradient 

with both train speeds and departure times variable, can be summarized as 

[ rj—(ckd} forie|l rij.} 

~{c{s,d}-il-?jX^-c7{di -ei;al)(7{di foiie{n^+l,...,2n^} 



(8-22) 



['dd,^ 

where reference is implicitly made to Equations 8-7, 8-12, and 8-15, 
Including the Economic Costs in the Gradient Search 

The types of costs incurred by raihroads for late deliveries were discussed above, 
and there was provided a differentiable approximation to the function of late costs 
expressed as a fimction of time. ' By using such an approximation, which is everywhere 
differentiable, the avoidance of late costs can be incorporated into the gradient search 
process. Arrival times are affected by both train speeds and train departure times, 
although either speed, departure time, or both may be varied during the search. 

The form of the late cost approximation function was given by (see Equation 7-:3) 

+ m,{u, -ti) for «, > 
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where 

Ui = the actual arrival time of the train, 
ti ~ the time at which late penalties begin to accrue, 
hi = the size of the step penalty (in k$), 
Mi = the rate of the linear portion of the penalty (in k$/hr.), and 
tc = the transition point where the cost function changes from a sigmoid to 
a line segment. 

Shortening this cost fimction to the form Ai(u^ for the nonce, and defining m = (wj ,...,m„ J, 
we may express a cost function which accounts for the arrival times of all trains in the 
form 

4«)=i:4(«,) (8-24). 

1=1 

But we also have the relationship 

u, = d, + - (8-25), 
so we may consider an alternative representation of Equation 8-24 as 

A(s,d)=ZAM (8-26). 

This latter form of the cost is appropriate to our search process, since that process 
is based on varying tiie components of the vector s and d . 

Now to incorpoiBte late arrival costs into fhe search, we extend the cost function 
of Equation 8-18 to the form 

D(s,d)=?],c{s,d)+T},E(^)+n,Af^,d) ' (8-27), 

where the 77. are weighting factors satisfying 

771+772+773=1 (8-28). 
The choices for these weights must be determined by experiment, and in one 
embodiment of the present invention it is possible to vary them iteratively as the search 
35 



wo 01/49550 



PCT/USOl/00080 



progresses. Individual users of the present invention may assign these weights as 
determined by the characteristics of the corridor and the costs imposed to the raibroad for 



weights take on values as determined in conjunction with the discussion of Equation (7-8) 
above. 

Gradient Search to Optimize Schedule Feasibility, Early Departures, and Economic 
Costs by Varying Only the Train Departure Times 

A search using the late arrival cost function of equation 8-27 may involve 

variation of only the departure times d , in which case the gradient by which the search is 
directed is of a form analogous to that shown in Equation 8-19. We may then express a 
component of the gradient vector in the form 



Borrowmg from Equations 8-6 and 8-20, we expand Equation 8-29 to the form 



the various effects built into the search algorithm. In the preferred embodiment, these 



(8-29). 



.(D{d4=r]\ 











(8-30), 
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which, with the help of Equation 8-7, provides aa exphcit representation of the 
components of the gradient of when only the train departure times are varied. 

5 

Gradient Search to Optimize Schedule Feasibility, Early Departures, and Economic 
Costs by Varying Only the Train Speeds 

If the dq)arture times of trains are held constant, and speeds are varied, then the 
10 gradient used to alter the speed vector s - (s, ,...,5,,^ ) during the search is of the form 

= 77,Vc(?,rf)+773v4?,^) (8-31), 

Where si^] is independent of the train speed. We therefore can obtain the k-th 

15 component of this gradient as 



= 77 A(c(j,4+%. 



(8-32), 



where the first term on the right side of Equation 8-32 can be expressed in a completely 
20 exphcit form by reference back to Equation 8-12. 

Gradient Search to Opttmize Schedule Feasibilily, Early Departures, and Economic 
Costs by Varying Both Train Departure Times and Train Speeds 

25 In this case, both d and s are variables in the complete cost function of Equation 

8-25, so the gradient takes the form 

V^''\D{^,d}=^Tj,V^^\ci^,d}+n^V^\E{^^^ (8-33). 
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Again regarding the gradient in vector form, the components of the first term of 
the sum on tiie left of Equation 8-33 are readily obtainable with the help of Equations 8- 
14, as explicitly represented witii the help of Equations 8-7, 8-11, and 8-12. The 
components of the second term can be obtained using Equation 8-20, and the components 
5 of the third term are obtained using Equations 8-30 and 8-32. 

Including Maximum Speed Limitation Effects in the Gradient Search 

A component of the cost function that would rise sharply in value as the speed s,- 
of a train Tt became close to Has maxunum speed sf^^ specified for the train was 
1 0 developed above. That component had the formulation (see Equation 7-9) 

= (8-34), 

and occurred as a weighted tem of the cost fimction, i.e., 

g(?,^)= J7,c(?,rf)+ 772^(?;rf)+ V,V^) (8-35) 

15 where the sum of the weights is chosen to be 1 in the preferred embodiment. Since 
variation of speed is independent of the departure times of trains, we have that 

V%(J) = 0 (8-36), 
so constraining the search by maximum train speeds does not effect the components of 
the gradient obtained as partial derivatives with respect to the departure times. 
■ 20 Relative to the gradient terms obtained as partial derivatives with respect to train speeds, 
we have 

= rjy^'^c[lsyn,V^Ui4^^^^ (8-37), 
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(8-38), 



where the exphcit form of the derivative in Equation 8-38 is from Equation A-2 of the 
Appendix, 

Expression of the Full Gradient 

For the sake of completeness, the complete expressions of these components of 
Equation 8-38 are provided below. First, let 



and select the weighting factors ?7i.'72>'73.'74 satisfying 

Note the indexing of the vector D places the partial derivatives with respect to s,, first 
and then the partial derivatives with respect to d,, second, but there are rij. values of each 
index. 




forA:e{l,...,«y} 
forite{«r+l 2«j.} 



(8-39) 






(8-40), 



where 
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¥oT ke{nj.+l,...,2nj.}, 
5 A=^.~ I {t[<^(y.-bj'tp~cr{y,-bj))~a{y,-ajjp~a{y,-aj))+ 

o{w-y,)(j3M^-y.))-<^(y.-e)iP~cr(y,-e))^^^ (8-42) 

-er>cclWi -m)) +73]^ for « > r 
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Illustration of the Gradient Search Process 

An example with twelve trains, six in each direction, running in a 150 mile corridor 
over an eight hour time window is provided below. 
5 Figure 13 shows the striag graph for the initial unprocessed schedule (i.e., train 

departure times were chosen without regard to feasibihty), and Table 1 below shows the 
information concerning each tram. There are twelve trains on the corridor and tiie time 
jframe of interest is eight hours (12:00 to 20:00). The colunms of the table indicate: 

(1) the train identification number (shown in the string graph as an integer at the center of 
10 each associated string) 

(2) direction of travel (Direction), 

(3) earliest acceptable departure time (Min. Departure), 

(4) actual departure time (Act, Dqparture), 

(5) latest arrival time before penalty is incurred (Max. Arrival) 
15 (6) initial speed (Speed), 

(7) train length (Length), 

(8) initial penalty incurred for being late (Penalty Step), 

(9) per hour penalty for each hour late (Penalty Slope), 

(10) maximum permitted speed (Max. Speed). 

20 
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The gradient search as discussed above was initiated with departure times being 
varied and train speeds held constant, and wifli the cost function including the penalties 
for early departure and economic penalties (i.e., late arrival). The resulting string graph is 
shown in Figure 14. 

5 Comparing Figure 14 with Figure 13, it can be seen that of the initial 31 points of 

intersection of train trajectories, in Figure 13, nine were close to feasible, where we will 
define "close" rather arbitrarily in terms of the intersection dots at least touching a siding 
20. Thus 23 intersection points were not close to feasible. In the final version of Figure 
14 some intersection points have disappeared, primarily because trains four and five have 

10 joined together in a convoy (identified in Figure 14 by the number 5 on the coincident 
strings), and some trains have been pushed off the string graph. In Figure 14, there are 
only two intersection points not meeting the definition of close. 

Table 2 below shows the final schedule, which resembles the original schedule 
except for the actual departure times of trains. Note that all trams require 7.5 hours firom 

15 actual departure until arrival at the destination, so only brain six is late, but train six is in 
fact only four minutes late. 
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5 Improving the Gradient Search Result by Speed Adjustments 

la this embodiment the gradient search result is modified by adjusting train speeds 
between sidings to achieve better siding meets. The gradient search process brought train 
intersections near, but may not have brought them exactly to the center points of sidings. 

10 This embodiment includes a technique for accounting for actual siding delays by 
changing intersiding speeds of trains as necessary to preserve the positions of intersection 
points at sidings. Ib order to provide a standard basis for that process, in this , 
embodiment we will first adjust the results of the gradient search so that the intersection 
pomts of train trajectories have y-coordinates precisely at the centerpotnts of sidings. The 

15 intersection points must be moved in order of increasing time coordinate, to assure that 
all prior intersection points have aheady been appropriately adjusted. 

To center intersection pomts at sidings and side specific trains, the train speeds of 
the trains involved must be modified somewhat. Of course, modifying a train's speed at 
any point could affect its trajectory downline, which would move the positions of its 

20 future meets with other trains. This is avoided by requuing that the centered intersection 
points remaui fixed, and tliat train speeds be varied as necessary to meet that requhement. 
More specifically, the tram that will not be sided at a given hitersection point will be 
constrained to pass through the centered intersection point, and the train that will be sided 
will undergo speed adjustments as needed to arrive and side before the opposed train is 

25 within an interfering (i.e., minimum stopping distance) of the siding train. 

The intersection points are processed in of increasing time order, so that all 
downline adjustments of trajectories may account for earUer modifications. As each 
intersection point is processed, the decision as to which train to side may depend on 
various criteria, which can be estabUshed as special rules auxiUary to the overall 

30 algoritimi. For example, if only one of the tsvo trains is too long for the siding, then the 
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other train must be sided. Another special case would be invoked for a train which could 
not restart if it sided on an upgrade of the corridor (that is, it could not generate sufficient 
tractive effort to move uphill). 

If there are no special circumstances dictating that one of the two trains should 
side, then the criteria for deciding the train to side is that of train speed; in effect, siding a 
train requires that it arrive "early" at tlie siding, relative to the centered intersection point, 
so that it can slow down and pull into the siding without interference from the opposed 
train. Arriving early implies that the train must obtain a speed greater than that which 
was nominally assigned by the gradient-search process of the present invention, and there 
is of course, some practical upper limit on train speed, as will be discussed below. The 
siding decision must be made based on which of the two trains will be driven less far 
toward its upper limit, given that it must be sided. Once the decision is made, the speed 
and arrival times of both brains are fitted to the actiial requirement of siding the brain. 

Figure 15 shows such a situation, where the intersection point (pcjj, ytj) of trains 21- 
and Tj is to be moved to the center of siding Sh (designated by point (xy, (a/, + b^) /2)), 
given that tlie immediately prior uitersection pomts affecting trains T, and Tj have ahready 
been adjusted. Clearly the speeds needed for train Tt (from Sh-i to S],) and for Tj (from 
jS/i+y to Sh) are given by 



(10-3) 



c-^lZ£l forh=l, 
* 2 



and trains li, Tp are the trains representing the immediately previous meets of with trains 
Ti and Tj, respectively. 
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There is also the case where tiaere is no prior point of intersection, i.e., where the 
intersection point {xij, yij) is the first intersection point for either or both of 7} or 2}, as 
shown in Figure 16. In this case, the speed needed to assure intersection at the siding 
center is given by 

5 

s - _Sl__ for Ti eastbonnd (10-4), 

and 

s„ = for Ti westbound (10-5). 

Xy-d. 

Figure 17 illustrates the result of centaing all meets for the gradient search results 
10 shown in Figure 14, on sidings 181 through 188 by adjusting train speeds between 
sidings. In effect, rather minor speed adjustments are usually sufficient to center all 
meets. 

Resolving Siding Conflicts 

15 

There is one possible undesirable side-effect that may arise when centering meets 
or passes, as illustrated in Figures ISA and 18B. After executing the gradient search 
process, the initial intersection points are shown in Figure 18 A. Train Tj intersects train 
at point 180, train Tj intersects train T4 at point 181 and train 7; intersects train at 
20 point 1 82. The result of centering all the meets represented by points 1 80, 1 8 1 , and 1 82, 
by speed adjustments as discussed above, is shown in Figure 18B. Train Tj is sided on 
siding Sn+i at point 183 because of its meet with train T3 , and train is sided at siding 
Sn+i at point 184 because of its meet with train . 

The difficulty created is that trains T] and T4 must both be sided on the same 
25 siding S„+i, although they are traveling in opposite directions, because one train is waiting 
on the siding that the other train must occupy before the former traia puUs out. This 
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cannot be accomplished, so the result of centering all of Ihe meets will in a case like this 
be an infeasible schedule. We will denote these artifacts as siding conflicts. 

The meet-centering process can produce two types of siding conflicts, as shown in 
Figure 19A and 19B. Figure 19A repeats the siding problem illustrated in Figure 18B. 

5 Figure 1 9B illustrates another siding conflict situation, but as in Figures 1 8B and 19A, the 
problem is again that two trains traveling in opposite directions must be sided on the 
same siding. Trains and intersect at point 194, with the former sided, while trains 
Ti and r4 intersect at point 196, with the former sided. Both of the siding conflict types 
shown in Figures 19A and !9B can be resolved by moving the meet of the conflicting 

10 trains to an adjaceaat sidmg, as shown in Figures 20A and 20B. 

Figure 20 A represents a siding conflict identical to Figure 19 A. The conflict at 
the meet point 200 is resolved by moving it upward to point 201 in Figure 20B. This is 
accomplished by accelerating or decelerating the necessary trains between adjacent 
sidings. Similarly, the siding conflict of Figure 19B can be resolved by moving it 

15 downward. 

This process of resolution as illustrated by Figure 20B (that is, the upward and 
downward movement of meets to resolve siding conflicts) will work if the trains in 
conflict have at most one meet at the siding to which their meet is moved, but will not 
work if both trains have meets at flie siding to which their meet is moved, as shown in 
20 Figures 21A and 21B. Li that case, the resolution of the original siding conflict at point 
210 in Figure 21 A by moving it to point 211 in Figure 21B, simply creates yet another 
siding conflict. 

However, there is an inductive way to resolve all siding conflicts which might 
occur from the meet centering process: if we call the siding conflict of Figures 18B and 
25 19A an upward-resolvable conflict, and the siding conflict of Figure 19B a downward- 
resolvable conflict, it follows that any siding conflict occurring on siding Si is in fact 
resolvable, because the conflict point may be pushed to the end of the corridor, where 
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any meets with the two trains involved can be avoided by slightly modifying the 

departure/arrival times of the involved trains as necessary. 

This is illustrated in Figures 22A-D, where the illustrations on the right provide 

resolutions of the siding conflicts on the left. The intersection at point 220 in Figure 22A 
5 is moved to point 221 in Figure 22B, by reducing the speed of train . hi Figure 22C, 

tiae siding conflict at point 224 is removed by moving the intersection of point of trains Tj 

and to point 225. Now feasible sidings caa occur at intersection points 225 and 226. 
Now we may proceed by induction to show that all siding conflicts are resolvable, 

with the basis beiug provided by the techniques demonstrated in Figure 22, and with the 
10 inductive assumption being that all siding conflicts occurring on siding Sn-u for ntl, 

can be resolved by pushing the conflict point to the end of the corridor. 

Figures 23A flurou^ 23E illustrate a downwaM-resolvable siding conflict on 

siding Sn, and it is shown that for all possible variations of that conflict, it can be resolved 

to a. situation where, at worst, it results in a new siding conflict resulting on siding Sn-i. 
15 By our inductive assumption, any such induced siding conflicts can be resolved. Figure 

23A shows the original meet situation. Figure 23B (case 1) shows the resolution if trains 
and have no meets at points d and f. Figure 23C (case 2a) shows the resolution 

when train S^has a meet at point d, train has no meet at poiut f, and train is not 

sided. Figure 23D (case 2b) ) shows the resolution when train r,has a meet at point d, 
20 train has no meet at point % and train is sided. The resolution of case 3 (not 

shown) where train 7; has no meet at pomt d, train has a meet at point f, is identical 

to case 2a and 2b. Finally, case 4 is illustrated in Figure 23E where trains and both 

have meets on siding S„_^ . 

Figure 24 shows a similar demonstration for an upward-resolvable siding conflict 
25 on S„, except the illustration is limited to a worst case, with it being evident that cases 

with fewer constraining meets are also resolvable to, at worst, siding conflicts on S„-j. 
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Figure 24A illustrates the original meet situation with the modification accomplished by 

moving the meet at point c to point g, as illustrated in Figure 24B. 

We conclude, finally, that although the meet centering process can produce 

infeasible string graphs because of the siding conflicts, all such siding conflicts can be 
5 resolved to feasible situations which do not include siding conflicts. When moving a 

meet point from one siding to the next lower one, there will usually be some horizontal 

latitude as to where to place it, and so to some degree, train speed limits can be favored. 

Note, however, that the resolution of these conflicts may resuU in some occasions where 

trains must travel at unrealizable speeds. This will be dealt with by iutroducing a new 
10 gradient optimization process in another embodiment of the present invention below. 

Accounting for Siding Time 

As described to this point, the invention permits an mitial schedule of trains on 
the corridor, arranged without regard for meets and passes, to be moved toward a 
1 5 schedule which minimizes or eliminates meets or passes occuiring at infeasible locations, 
i.e., not at sidings. 

After the processes of improving the gradient search results by speed adjustments 
and resolving siding conflicts have been apphed, as discussed above, to the original 
gradient seai-ch result, there has been created a string graph m which each train trajectory 

20 is depicted as a sequence of straight hne segments, constrained to meet other train 
trajectories at the centerpoints of sidings. The string graph, adjusted after the gradient 
search as necessary to move all meets to centopoints of sidings, will be called the 
incomplete string graph. 

The gradient search and the speed adjustments produce a meet of two trains at a 

25 siding, but in one embodiment it does not actually account for the need for one train to 
side, or for the fact that the train has a length. To actually side one train, it must arrive at 
the siding far enough in advance of the other train to completely pull into tiie siding, and 
it must delay its departure until the other train is clear of the siding. 



50 



wo 01/49550 



PCTAJSOl/00080 



Figure 25 illustrates this problem. To this point, a train trajectory has been 
approximated as a single unbroken line segment (as in Figure 2), in actuality, it will take 
the form of a broken line segment if tiie corresponding train must be sided. In Figure 25, 
the trains T2 and T3 must side, so the corresponding trajectories L2 and L3 reflect the 
required siding time with horizontal line segments 250 and 252 inserted into the 
trajectories. The minimum length of the horizontal segment is determuied by the length 
and speed of the opposed tram. Therefore the level of resolution into train trajectory 
planning must be improved m this embodiment to obtam an unplementable train 
schedule, based on the results of the gradient search. It is necessary to develop the 
mathematics of siding trains, given that an initial schedule has been obtained using the 
gradient search process, above. 

Defining the Train Trajectory Vector 

Implicit in the purely geometric format described to this pomt, are the numerical 
quantities needed to define the train trajectory vector of Equation 10-1. Specifically, for 
tram Tt, the value of bw {Ti eastbound) or of ,,^,+1 {Ti westbound) must be equal to the 
departure time di of the train, which was determined by the gradient search process, with 
possible modification by the resolution of sidmg conflicts. Now for an eastbound train, 
assume that the first meet with another train occurs with traui Tj at siding 5a, ^ > 0, thus 
we specifically know that Tt must be at pomt {xy, a) on the string gr^h, as shown in 
Figure 26. Then the speed Sih of train Tj, firom its origin, must be 

and it follows that bik , for k = 1, h, and eik, for k = l,...h-l, may be determined as 
follows. 

**=*/o+— » for A:=1,...,A (lO-U 

and 
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%=^,o+" forA:=l,...,M (10-8). 

We can now proceed on the next line segment {i.e., from meet to meet) defining 
the trajectory of 2} to obtain a speed, determined by the intersections of Ti with other 
trains, from which we can detamiae the times of arrival of 2} at all intermediate siding 
5 edges, thereby filling ia all of the data required for the train trajectory vector of 2} except 
the siding decision values 5,/,. Siding decisions have not yet been considered, so these 
values vidll be defiined later. 

It should be clear that an analogous process can be defined for westbound trains, 
so we have inductively defined all traia trajectory vectors using the incomplete string 
10 graph. 

Extending the Definition of flie Train Trajectory 

The definition of train trajectories as equations relating distance along the corridor 
to time, as given by Equation 3-3, does not accommodate the siding time and siding 
15 decisions required for some trains. Instead, it provided a characterization of trajectories 
as straight line segments, for the purpose of miiiiinizing the computations needed for the 
gradient search process. In order to generalize the trajectory, in this embodiment the 
simple definition of a trajectory will be modified by adding parameters accounting for 
tram delays at sidings. 

20 Where a corridor has ns sidings, we begin by defining the train trajectory vector, 

and for notational convenience, we will designate the west end of the corridor as siding 
So, and the east end of the corridor as siding S„^^i, with the recognition that these 
"sidings" have a length of zero. Given this convention, define the train trajectory vector 
for train Ti as 

25 M,-(d,AoA 5,i....,5,J (10-1), 

where 

= the direction of train Ti (akeady defined in Equation 3-4), 
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bih = the time at which train Ti reaches siding Sh (h = 0,.. .,ns+i), 
en, = the time at which train Ti departs the siding Sh{h = l,.. .,ns). 




ri if r, is sided on 5a 

1^0 if Ti is not sided on 



The times at which a train reaches or departs from a siding will be the time at 
5 which the head of the train reaches the upstream or downstream ("downstream" or 
"upstream" is defined relative to the direction of the train) end of the siding, respectively. 
Also for consistency, since siding Si has endpoints a,- and bi, as measured from the west 
end of the corridor, let bo denote Ihe beghining of the corridor, and a^^+i denote the end 
of the corridor. 

10 

Detailing the Siding Process 

Figure 14 demonstrated the result of the gradient search process, and demonstrates 
that the search process has the capability to adjust departure times so that train trajectories 
intersect at sidings. The gradient seai-ch process cannot usually perfectly align all meets 
15 at sidings, and so the meet-centering process was also described above. Once we have in 
fact placed all meets at sidmgs, we might use train speed adjustments to interpret the 
resulting string graph as showmg that the engines of trains pass exactly at the centerpoints 
of sidings. 

Now the focus will be on a technique by which a string graph such as that m 
20 Figure 14, with meets centered on sidings, as discussed above, can be modified further to 
provide a full, feasible string graph schedule with trains sided as necessary. We will 
assume ttiat we begm with all train meets centered at sidings, and all possible siding 
conflicts resolved as necessary. The process vsdll be inductive: we will begin by ordering 
the collection of all intersection points on the incomplete string grs^h according to the 
25 time of intersection, and we will proceed to modify them, in time order, so that each 
intersection point reflects a feasible siding arrangement. 
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Figure 27 illustrates the technique to be applied, as it is applied to intersection 
point y24- It is assumed that all intersection points of the string graph prior in time to y24 
have abready been modified by this process, so that flie required time and speed data 
concerning trains T2 and T4, prior to point y24 are in fact valid. Of the two trajectories 

5 passing though y24, we will choose to side train T4, and the modification of trajectory for 
T4 is indicated by the dashed sequence of Hne segments. Effectively, we require that T4 
operate at a higher speed from the last intersection pomt on the trajectory (relative to the 
incomplete string graph) in order to arrive at siding Sn so that the last car of T4 actually 
enters the siding before the engine of train T2 arrives at the west end of siding 

10 Figures 28 and 29 represent possible meet/pass situations between trains. There are 

four basic cases, as follows: 

(1) an eastbound train sides for a westbound train, 

(2) a westbound train sides for an eastbound train, 

(3) an eastbound train sides for a passing eastbound train, 
15 (4) a westbound train sides for a passing westbound train. 

There are also four variants on each case (for a total of 16 cases), depending on 
whether either or both of the trains mvolved were sided at the previous intersection point 
on their trajectories. This has significance because a train leaving a siding Avill have a 
lower initial speed (the pullout speed firom the siding) across an intersiding segment than 
20 a train which has not been sided. 

Essential parameters for the process will be defined in conjunction with Figures 28 
and 29. Relative \o any train 2/, let 
Aih - the arrival time of the last car of T,- at the upstream end of the sidmg Sk, 
Dih = the time at which the head of T; arrives at the downstream edge of siding Sb, 
25 tih = the time at which a train Ti not sided at Sh passes the midpoint of the siding Sh, 

Vk = the pullin/pullout speed of any tram for siding Sh, 

p(i,h) = the siding at which Ti had the most recent meet before the present meet at 5/, . 
fi(v) - the minimmn stoppmg time for train Tt at speed v. The approximation used for 

this function is explained in Appendix B. 
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Although the following notation is not new, we review it here for convenience: let 
ah = the coordinate of the western end of sidiag iS/„, 
bh = the coordinate of the eastern end of siding Si,, 
Mi = the length of traia 3}, 

L = the length of the corridor (with the western end being the origin), 
and finally, let 

c,^ = Ejilh. ^ the coordinate of the midpoint of siding Sh (1 0-9). 

Relative to the earlier descriptors of the train trajector vector for train Ti (Equation 

10-1), note that 

DiH^Bi, (10-10), 
and, for an unsided train passiug siding Sh, assumed to mamtaiQ constant speed across the 
extent of the siding, 

t., ^hlliL (10-11), 

the value of which was estabUshed by centering all meets at sidings. 

In the following derivations, the trains meeting at a siding Sh vwll be trains Ti and 
Tj, and 2/ will always be the traui to be sided. The apparent constraints that must be met 
for Ti to be sided (see Figures 28 and 29), are 

\^Dj, (10-12), 

and 

D,>Aj, (10-13). 
These two constraints are somewhat idealized, and both need modifications. 
First, it would be unsafe to apply inequality 10-12 Uterally, because if, for any reason, 
train 1/ were to stop short of being fiilly sided, then train 3} might hi fact be too close to 
stop in time to avoid a collision. Thus condition 10-12 should be rq)laced with 



55 



wo 01/49550 



PCT/USOl/00080 



where vj is the speed of 2} as it approaches siding iS*. 

Condition 10-13 also requires a modification, because it could be the case that 2} 
could actually clear the downstream end of the siding (relative to Tt) before T/ could get 
there, evoa if TJ pulled into the sidiag, continued moving at maximum siding speed, and 
5 arrived at the downstream end of the siding, hi that case, Doi is limited by die speed of 
Ti, not the position of Tj, and takes the minimum value 

j^.,^^4^_4- ^^-^;"^' (10-15), 

so the corrected version of Condition 10-14 is 
10 iP,=max|4,.A + ^'~^^^"^' | (10-16). 

Constraints 10-14 and 10-16 then provide the practical constraints by which meets 
and passes can be planned. 

The quantities in the mequalities are functions of the train speeds on the previous 

15 intersiding segments and of the departure times from the last siding: inductively, we 
assume that the departure times for both trains from their previous meets are known, and 
we must derive the speeds needed by both trains to anive at the siding Sh so that 
constraints 10-14 and 10-16 are met. The known quantities, for the trains T,- and Tj, at the 
beginning of the inductive step, are 

20 (1) tjh, the time at which 2} should be at the center of 5/, (Equation 10-1 1), 

0) Dj.po.it) for Tj. 

To satisfy constraints 10-14 and 10-16, we must detennine values fox Da,, Djh, An, 
and Ajf„ in terms of speeds, and then solve the constraint inequalities for the speeds 
25 required to meet the constraints. 

The speeds so obtained are for 1/ and Tj from their last meets to their common 
meet, and when we solve the constraint mequalities (subject to the siding choices made) 
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to obtain these ttain speeds, we will also determine the values corresponding to items (1) 
- (3) above for trains 1} and 2} at siding Sh, thereby completing the inductive process. 
The basis of the induction will be taken up later. 

The Inductive Step for the Unsided Train 

We first determine the speeds, from the requirements that the unsided train pass the 
center of siding Sh at time tju: 

c —b 

D, + Ji — EiMl for Tj eastbound, and not sided at S^^^j^^^ 





. + ^(^•■*> for Tj eastbound, and sided at 



for Tj westbound and not sided at S^^j^^^ 



— i for Tj westbound and sided at S,, 



(10-17). 



Note that thore are also two special cases of Equations 10-17, namely 



for Tj eastbound. 



(10-18) 



for Tj westbound 



and for notational consistency, we define 



and 



(10-19). 



= L 



57 



wo 01/49550 



PCTAJSOl/00080 



Additionally, the validity of equation 10-17 requires that the distance between 
sidings Sh and Spffj,j exceed the length jy, of train Tj. 
From Equations 10-17, we may solve for the speeds required: 

= ^ ^'•"^pu.h) ^ Jqj. j.. eastbound and not sided at Sp(i,h) (10-20); 



■i— , for Tj eastbound and sided at Sp^j,) (10-21); 



s.^ ^ Jim — i- , for Tj westbound and not sided at Sp(f,h) (10-22); 

s = ^PU.i>)-^u~Mj westbound and sided at Sp(j,H) (10-23). 

■"' Mj 

^p{m 

Now that the speed for the unsided train is determined, we may solve for Djh and 
Ajh, as follows: 

r b -c 

tji, +— — - for Tj eastbound 
D„=\ . (10-24); 



for T, westbound 



— "'•'^^J fot3\ eastbound 

f — * — L forr, westbound 

J" ^ ■' 



(10-25). 



For the unsided train, the determination of Djh in Equation 10-24 completes the 
inductive step of the siding algorithm. Note that if there are sidings between Sh and 
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Spo\h), then the times of arrival and departure from those sidings is implicit in the speeds 
calculated in Equations 10-20 through 10-23. For k an index of such a siding, we have 
tiie following results: 



for Tj eastbound and sided at Spg^^^ 
Tj westbound and not sided at S^g ,,) 



3^^^^.^^ + J£l. + ^pu.") ^* for Tj westbound and sided at Sj,g 



We may then write 



5 4. h — 2i. for Ti eastbound 



- for r, westbound 



(10-26). 



(10-27). 



The Inductive Step for the Sided Train 

Half of the inductive step for the sided train Ti is already complete, in that we can 
set the value D,/, to be any value satisfying condition 10-16, although we would normally 
set that value to be as small as possible. However, we must also determine the speed 
required by I}, from the previous siding Sp^j,) where Ti had a meet to S;,, that will satisfy 
Condition 10-14. There are four cases, based on whether T, is eastbound or westbound, 
and did or did not side at /,;. We express Aih for each of these cases, and then use 
Condition 10-14 to determine a minimum speed for Ti, 
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+ ^pW'^ + Ml for 2] eastbound and not sided at S^^^ 

+ J£j_ + ^H-b,(i,,)-M, ^ gjj. ^ eastbound and sided at S^^, 
+ "piW^'' +Ml for 2; westbound and not sided at ^^y^j^ 

. J£l. + ^(W-^ft--^. + jA^L for T. westbound and sided at S,,, « 



(10-28). 

Since the value of Dji, has been determined in the previous section. Equation 10- 
28 and Condition 10-12 lead to inequalities for the speed sut or Si^-i of 7/, as follows. 



- , for Ti eastbound and not sided at Sp(t,h) (10-29), 



- , for Ti eastbound and sided at Sp(i,h) (10-30), 



.Ml 



, for Ti westbound and not sided at Sp(i,i0 



(10-31), 



- , for Ti westbound and sided at Sp(i,h) (10-32), 



, 10 



[ for 2} eastbound 
. r, westbound 



(10-33). 



All of the quantities on the right sides of inequalities 10-29 through 10-32 are 
known, so the speed Sik or Si^k-i^ox train Ti is determined, and the haductive step is 
complete. If a siding Sk is intermediate to Sh and Sp(i.h), then Equations 25 and 26 
establish the values of e/t and bit. 
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Establishing an inductive basis for the above depends only on the observation that 
the very first meet for any train Ti or 2} is preceded by the entry into the corridor firom the 
east or west end. All of the computations then required to arrive and meet at siding Sh, 
subject to the constraints, are based on the original departure time of the relevant train, to 

which Dpf/,/,; orDp(fji) is set equal, as the case maybe. 

Finally, the inductive process defined above determines speeds, and the times of 
arrival and departure for each train at each siding, based on meets at the sidings. Once a 
train has encountered its last meet, the final speed is adjusted to assure that it arrives at 
the end of the corridor as scheduled. If train T, has its last meet at siding Sh, then the 
speeds between all subsequent sidings which are required to exit the corridor as 
scheduled are given by 

L-b, 



L-b^-Mj 



- for Tt eastbound and not sided at 



-Ml 



for Tf eastbound and sided at 



(10-34), 



0;, -M, 



for Tj westbound and not sided at 



.Ml 



for Tj westbound and sided at 



(10-35). 



Figure 30 displays a final and complete string graph which was has been adjusted 
for centered meets, and then for train sidings. 

Figure 31 is a flow chart implementing one of the algorithms of the present 

invention. The flow chart of Figure 31 can be processed on any special purpose or 
general purpose computer. The software code necessary to implement the Figure 31 flow 
chart can be written by anyone who is skilled in the art of preparing software code, given 
the information in Figure 31 and the description of the invention provided hereui. 
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Processing begins at a step 310 where tiie initial conditions are established. That 
is, there is assumed an initial vector)?, , which identifies either the initial train speed or 
the initial departure times of the trains on the corridor, or both. The vector y„ is used to 
calculate the intersection points at a step 312 and then the value of the localizer function 

5 for each calculated intersection point is determined at a step 314. At a step 316, the 
localizer ftmction values are summed to create a schedule feasibility cost function with an 
argument y„. As discussed above, there are many different cost fimction types 
associated with different embodiments of the present invention. For instance, Equation 8- 
17 identifies two cost functions. The cost ftmction of schedule feasibility (Q and a cost 

10 function associated with early departure effects {E). The economic cost function is 
defined in Equation 8-26 and the maximum speed cost function is defined in Equation 7- 
9. Depending upon the embodiment of the present invention, one or more of these cost 
functions will be used to create the cost function at the step 316. 

At a step 318, the gradient of the cost fimction at y„ is calculated. At a step 320, a 

1 5 new argument for the cost function is created. This argument is referred to as y„^^ and is 
calculated using the gradient value fix)m step 318 and a predetermined step size. This 
step size is based on the gradient value and must be determined in each situation so as to 
converge toward the fimction minimum. Reference is made to the four step process 
outlined at Equation 8-3. At a step 322, the magnitude of the difference between the cost 

20 function at y„ and is calculated. At a decision step 324, the results from the step 
322 are compared to a flireshold. If the threshold is not exceeded, then the cost fimction 
minimum has been located and a schedule for the corridor is produced. This is illustrated 
diagrammatically at a step 325, If the threshold is exceeded, then further calculations can 
be performed to find the cost function roinimum. At this point, processing moves to a 

25 step 326 where the previous value of 7„ is now set equal to the value of 7,,+, and 
processing returns to the step 312 where the intersection points are again calculated. 
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Processing &en continues throu^ the steps 314, 316, 318, 320, and 322, followed by 
decision step 324 where the magnitude is again compared to the threshold value. 

As discussed above, there are additional refinements that can be made from the 
5 schedule produced at the step 325. These refinements represent additional embodiments 
of the invention and are discussed in detail above. In flow chart form, they are presented 
in Figure 32. hi lieu of processing proceeding to the step 325 in Figure 31 when the 
threshold value is not exceeded, processing can instead continue to a step 340 illustrated 
in Figure 32. Here, adjustments are made to intersiding train speeds so that the 
10 intersections will occur precisely at the sidings. This embodiment is discussed in 
conjunction with Figures 15, 16 and 17. In another embodiment, siding conflicts can be 
resolved at a step 342. This embodiment is discussed in conjunction with Figures 18—24 
above. The matter of accounting for the time the trains are sided is represented by a 
processing step 344. This embodiment is discussed above in conjunction with Figures 
15 25—30. Finally, incorporation of these additional embodiments provide for the 
generation of another train schedule for the rail corridor, as illustrated at a step 346. 
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Appendix A - Properties of the Sigmoid and Localizer Functions. 

Lemma 1 : a{x;a,fi) =fi-a{- xia,fi) (Al) 
5 Proof: 

Q.E.D. 

Lemma 2: -^(a{ax + b;a,jS))=^<T{ax + b;a,p'){fi -cT{ax+b;a>p)) (A2) 
ox p 

Proof: 

10 ^(cxica + b;a,p)) = -^{^{l + e-"<«^^>)r') = l)(l + e'^^-^'^Yi- aa)e-"^-^'^ 

= ^a{ax + b;a,/3)a{- {ax + b^a,p)=y-a{ax + b;a,jBjj3 - cr{ca + b;a,fi)) 

Q.E.D. 

15 The following result follows from the derivation of the localizer function in the body of 
this document, and from application of Lemma 2. 

Lemma 3: Let (-«3,w),(ai,&i),...(a^,&^),(e,oo) represent mutually disjoint intervals, 
with - 00 < flj < 6i < a^... < Ojv < < °° • Th®^ L{x\a,P) defined to be low if and only if 
X is in or near one of the intervals (-00, w),(aj,&,),...(a^,6;y),(e,<») takes the form 

20 Lix;a,fi) = y8 - 'f^<r(x - b,;a,fi)- a{x - ;(?)]- cr(w - x;a,p)-<T{x - e',a,l5), 
and 
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-a\x-aj;a,p\p-a\x~a;\a,l. 



A(i(^;a,,3)) = ^{|[o-(x-Z,,;a.^)0s-a(x-6,;a.;S))- 

+ a{w-x;a,fi){p - a{w - x;a,fi))~ a{x - e;a,pifi - a{x - e;a,/?))} (A3) 
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Appendix B - A Train Stopping Time Approximation 

The basic fonnula for acceleration/deceleration of a body is 

F = MA (Bl), 

5 where 

F = the braking force applied, 

M = the mass of the body, 

^ = the acceleration of the body. 

A train has brakes on every car, and each car has mass, so we will assume that the total 
10 maximum braking force and mass are proportional to the length of the train. Therefore 
Equation Bl may be written as 

A=k (B2). 
i,e., the deceleration available at maximum braking is (approximately) independent of ttie 
train's length or mass. 

15 To evaluate k, we assume that a train moving 50 mph could stop m 1 mile, therofer its 

average speed during (linear) deceleration would be 25 mph, and the time required to 

reach a full stop would be 

(1 mi. /25 mph)(60 min/hr) = 2.4 minutes. 

Thus the equation relating train speed v to stopping ismef(y) takes the form 
JO f(v)=v/A-=v/k (B3). 

and 2.4 = 50/k, 

tiierefore 

= 50/2.4 = 20.83 (mph/min). 

The final form is then 

25 /(V;=v/20.83 (B4), 

where f(v) is in minutes, and v is in miles per hour. 
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Appendix C - Variable Glossary 

Ai(t) - the late penalty function assessed for train Ti (Equation 7-1) 

Aih - the time of arrival of the end of train Ti at the upstream edge of siding Sj, 

5 A^, d)- the cost fimction component forcing on-time arrivals 

ah - the distance ftom the west end of the corridor^ at which siding Sh begins 
Bih - the decision variable as to whether train T/ is sided at siding Sh 
b/i - the distance from the west end of the corridor at which siding Si ends (fl/ < bi) 
c{s,d)- the cost function forcing trajectory intersections at sidings 

10 Ch - the midpoint of siding 5/, 

Dji, - the time of departure of train Ti from the downstream edge of siding Sh 
d{ - the departure time of train Ti 

d — the vector of dimension of departure times for all trains 
E - the name of the point at the east end of the corridor 
15 -the cost function component preventing early train departures 

fify) -- the minimum stoppiag time of train Ti at from speed v 
g(?,^)~ the total schedulhig cost function (Equation 7-10) 
Hi - the length of siding Si 

hi - the step penalty cost incurred when train TJ arrives late 
10 I - the set of all intersections of train trajectories (even if not on the string graph) 
Z ~ the length of the corridor 

Li - the line on the string graph representing the trajectory of train Ti 
L(y) - the localizer fimction, with minima corresponding to each siding (Equation 5-5) 
i(y) ~ the balanced localizer function (Equation 5-7) 
i5 Lij(yi^ - the localizer modified so trains Ti and 2} won't meet where neither can side 
Mi ~ the length of train Ti 
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nti - the late penalty per time unit when train Ti arrives late 

ns - the number of sidings along the corridor 

TiT — the number of trains involved in the optimization 

p(i,h) - the siding prior to Sh at which tram Ti had a meet, 

5 Si - the designator for the i-th siding, travehng eastward on the corridor 

Si - the speed of train Ti 

sf^^ -maximum speed permitted for train 2} 

Sih - train T/'s speed between the downstream edges of sidings S% and 

s -the vector of dimension Kr of speeds for all train 

10 T{ - the designator for the i-fh train 

Ty ~ the set of aU sidings where at least one of trains 1} and 3} can side 

ti ~ the arrival time at which train Ti begins to incur late penalties 

tjh - the time at which train Tj, if not sided at S;,, reaches q. 

tij - the time coordinate associated with traj ectory intersection point yy 

1 5 V{s) - the cost function component which lunits train speeds 

Vh - the pulIin/puUout speed for trains at siding Si, 

W — the name of the point at the west end of flie corridor (zero on the distance axis) 
— the distance j&om the west end of the corridor at which trains 2} and 2} intersect 

y -- the vector of all traj ectory intersection points yy 

10 a — the sigmoid function parameter controlling steepness of rise (Equation 4-1) 

p ~ the horizontal asymptote of the sigmoid function (Equation 4-1) 

7jj —the weight applied to the feasibility component c(j, d) of the cost fxmction 

rj2 — the weight applied to the late arrival component , j) of the cost function 

r}3 — the weight appHed to the early departure component E[d) of the cost ftmction 

15 T]4 - the weight applied to the maximum speed component V{s} of the cost jflmction 

0i - a variable denotimg the direction of brain Ti, 0 if eastbound, 1 if westbound 

cr(x) ~ the sigmoid function Equation 4-1) 
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WHAT IS CI.AIMED IS: 

1. A mettiod for scheduling the movement of a plurality of trains operating on a 
rail corridor, wherein trains traversing the rail corridor may intersect on the same track, 
whereby each train has at least one travel parameter that can be varied, whereby the rail 

5 corridor includes at least one main line and a plurality of secondary tracks onto which a 
train may be moved to avoid an intersection with another train, said method comprising 
the steps of: 

(a) deriving a localizer function to represent the rail corridor, wherein said 
localizer function has a value in a first range between secondary tracks and has a value in 

[ 0 a second range in the vicidty of each secondary track; 

(b) selecting a value for at least one travel parameter for each of the plurality of 

trains; 

(c) finding the intersection points for the plurality of trains; 

(d) determining the value of said localizer function for each mtersection point; 

L 5 (e) combining said localizer function values to create a cost function value; and 

(f) changing one or more of the values selected in step (b) to jSnd the roiiiimum 
cost function value. 

2. The method of claim 1 whereui the travel parameter includes train speed. 

3 . The method of claim 2 wherein the speed of at least one train is varied between 
JO intersection points to ensure the intersection occurs at a secondary track. 

4. The method of claim 1 wherein the travel parameter includes the entry time of 
the train onto the rail corridor. 

5. The method of claim 1 wherein the travel parameter includes tram speed and 
the entry time of the train onto the rail corridor. 

}5 6. The method of claim 1 wherein the secondary track includes a passing siding, 

7. The method of claim 1 wherein the secondary track includes two parallel tracks 
with crossover switches therebetween. 

8. The method of claim 1 wherein the localizer function is derived by summing a 
plurality of sigmoid functions, wherein said sigmoid fimctions are translated and inverted 
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with respect to each other and the location of flie secondary tracks such that the localizer 
function has a value in the jgrst range between secondary tracks and has a value within the 
second range in the vicinity of each secondary track. 

9. The method of claim 1 including the step of forming a first vector having as its 
5 elements the values selected in step (b), and wherein step (f) includes: 

(fl) determining the gradient of the cost function at tile first vector; 

(f2) calculating a second vector by subtracting the product of a predetermined step 
size and the gradient calculated in step (fl), from the first vector; 

(f3) calculating the cost function value at the second vector; 
iO (f4) calculating the magnitude of the difference between the cost function value 

calculated in step (e) and the cost fimction value calculated in step (f3); 

(f5) if the result of step (f4) is greater than a predetermined threshold, returning to 
step (fl) wherein the second vector is now used in lieu of the first vector in the gradient 
calcxilating step (fl). 

15 10. The method of claim 1 wherein the localizer function of step (a) is modified 

to account for rail corridor end points by setting the localizer fimction value to a value 
within the first range beyond the end points of the rail corridor. 

11. The method of claim 1 wherein step (e) includes the process of su3nming said 
localizer fimction values to create the cost function value. 

JO 12. An apparatus for scheduling the movement of a plurahty of trains operating on 

a rail corridor, wherein trains traversing the rail conidor may intersect on the same track, 
whereby each train has at least one travel parameter that can be varied, whereby the rail 
corridor includes at least one main Une and a plurahty of secondary tracks onto which a 
train may be moved to avoid an intersection with another train, said apparatus 

15 comprising: 

means for deriving a localizer fimction to represent the rail corridor, wherein said 
localizer function has a value within a first range between secondary tracks and has a 
value within a second range in the vicinity of each secondary track; 
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means for selecting a value for at least one travel parameter for each of the 
plurality of trains; 

means for finding tih.e iatersection points for the plurality of trains; 

means for determining the value of said localizer function for each intersection 

5 point; 

means for combining said localizer fimction values to create a cost function value; 

and 

means for changing one or more of the selected values to find the minim um cost 
function value. 

10 13. The apparatus of claim 12 wherein the travel parameter includes train speed. 

14. The apparatus of claim 13 including means for adjusting the speed of at least 
one train between intersection points to ensure that the intersection occurs at a secondary 
track, 

15. The apparatus of claim 12 wherein the travel parameter includes the entry time 
15 of the train onto the rail corridor. 

16. The apparatus of claim 12 wherein the travel parameter includes train speed 

and the entry time of the train onto the rail corridor. 

17. The apparatus of claim 12 wherein the secondary track includes a passing 

siding. 

10 18. The apparatus of claim 12 wherein the secondary track includes two parallel 

tracks with crossover switches therebetween. 

19. The apparatus of claim 12 wherein the locaUzer function is derived by 
summing a plurality of sigmoid functions, wherein said sigmoid functions are disposed 
with respect to each other and the location of the secondary tracks such that the localizer 

15 function has a value within the first range between secondary tracks and has a value 
within the second range in the vicinity of each secondary track. 

20. The apparatus of claim 12 including: 

means for forming a first vector having as its elemaits the selected values, and 
further including: 
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means for detennitiiag the gradient of the cost fuaction at the first vector; 

means for calculating a second vector by subtracting the product of a 
predetermined step size and the gradient calculated, &om the first vector; 

means for calculating &e value of the cost function at the second vector; 
5 means for calculating the magnitude of the difference between the cost function 

value based on the first vector and the cost function value based on the second vector; 

means for determining if the magnitude of the difference is greater than a 
predetermined threshold; and 

means for recalculating the gradient when the magnitude of the difference is 
[0 greater than a predetermined threshold, wherein the second vector is used in heu of the 
first vector in the gradient recalculation. 

21. The apparatus of claioa 12 wherein the localizer function is modified to 
account for rail corridor end points by setting the localizer function value to a value 
within the first range beyond the end points of the rail corridor. 
[5 22. The apparatus of claim 12 wherein means for combining includes means for 

summing said localizer function values to create the cost function value. 

23. A method for scheduling the movement of a plurality of trains operating on a 
rail corridor, wherein trains traversing the rail corridor may intersect, whereby each train 
has at least one travel parameter that can be varied, whereby the rail conidor mcludes at 
10 least one main line and a pluraUty of secondary tracks onto which a train may be moved 
to avoid an intersection with another train, said method comprising the steps of: 

(a) selecting a value for at least one travel parameter for each of the plurality of 

trains; 

(b) finding tiie intersection pomts for the plurality of trains; 

25 (c) assigning a first value to each intersectiion point occurring at or near a siding 

and assigning a second value to each intersection point that is not near a siding; 

(d) combining said values firom step (c) to create a cost fimction value; and 
(e) changing one or more of the values selected in step (a) to find the minimijm of the 
cost function. 

72 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCTAJSOl/00080 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USOl/00080 



2/16 





(0.0) 

FIG. 5 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCTAJSOl/00080 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USO 1/00080 



4/16 

^80 




W S1 S2 E 

FIG. 8 




FIG. 9B 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USOl/00080 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCTAJSO 1/00080 



6/16 



FIG. 12 




E 

150 
140 
130 
120 
110 
100 

Z 80 

5 60 

50 
40 
30 
20 
10 
W 



STRING GRAPH 















< V 




- v 


XX ^ i 




/ X / \ ,12- \ 1 




\,Xx,xX, 1 



20 



20 



M2 13 14 

F/G. 13 



15 16 17 
TIME 



18 19 



20 



SUBSTITUTE SHEET (RULE 26) 



PCTAJSOl/00080 



E STRING GRAPH 




FIG, 14 TIME 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCTAJSOl/00080 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USOl/00080 




T2 'T3 

FIG. 19A 



FIG. 19B 




T1 T4 T3 T2 

FIG. 20 A 



Ti T4 T3 \ 

FIG. 20B 




FIG. 21 A 



FIG. 21B 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USO 1/00080 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USOl/00080 



11/16 




F/a 23 A 



3^ 




T| ^5 T2 

FIG. 23C 



V ■ V 



F/G. 23D 




^1 Ts' Te ^2 
FIG.23E 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



12/16 



PCTAJSOl/00080 




Figure 26: 



wo 01/49550 



PCT/USOl/00080 



13/16 



t>3- 
33- 



ai - 



W 

FIG. 27 





T2\ 


/ '•'3 T4 y 














\ / ' • 
\ / ^ • 

V24// 






— _\^// — j 


/ \ '/ \ 
/ \ ' / \ • 
/ \ '/ \ • 


/ \. > yi4 \ 




^4/ 


\ V2 



TIME- 



\t 


h+1 f 


SIDING Sh+1 ] 


Di,h+l' 
TRAIN 1 


\\ /V^^TRAIN Tj 




V SIDING Sh 






/Ajh D|h 






SIDING Sh-1 



TIME- 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USOl/00080 




SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USO 1/00080 



FIG, 31 



15/16 

( START ) 

ASSUME VECTOR OF INITIAL 
TRAIN SPEEDS AND/OR 
DEPARTURE TIMES 
Tn 



—310 



CALCULATE INTERSECTION 
POINTS 



-312 



DETERMINE VALUE OF 
LOCALIZER FUNCTION FOR 
EACH INTERSECTI ON POINT 



—314 



SUM LOCALIZER FUNCTION 
VALUES TO CREATE COST 1 
FUNCTIONAT % "^"^'^ 



I 



COMPUTE GRADIENT OF 
COST FUNCTION AT Tn 



-318 



COMPUTE NEW ARGUMENT OF 
COST FUNCTION Tn+l 
USING GRADIENT AND 
PREDETERMINED STEP SIZE 



I 



-320 



COMPUTE MAGNITUDE OF 
DIFFERENCE BETWEEN COST 
FUNCTIONS AT AND Tn+1 

i 

DOES MAGNITUDE EXCEED 
A THRESHOLD 
? 



-322 



324 
NO 



325 



PRODUCE TRAIN 
SCHEDULE 



SUBSTITUTE SHEET (RULE 26) 



wo 01/49550 



PCT/USOl/00080 



16/16 



( START ) 

ADJUST INTERSIDIN6 
TRAIN SPEEDS 



I 



RESOLVE SIDING 
CONFLICTS 



ACCOUNT FOR TIME 
TRAINS ARE SIDED 

i 



PRODUCE TRAIN 
SCHEDULE 



-340 



"342 



-344 



"346 



FIG. 32 



SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



Inte lal Application No 

PCT/US 01/00080 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 B61L27/00 G06F17/60 




Accordinq to Internationa) Patent Classification (IPC) or to Irotii naUonal classification and IPC 




B. FIELDS SEARCHED 


Minimum documentation searched (classlficalion syslem toliowed by classification symbols) 

IPC 7 B61L G06F 605D G06N G08G 


IJocumentation searched other than minimum documentation to the extent that such documents aie Included In the 


ieids searched 



Eleclronic data base consulted during the international search (name of data base and, where practical, 

EPO-Internal , WPI Data, PAJ 



C. DOCUMENTS CONSIDERED TO BE RELEVAWT 



Relevant to Claim No. 



US 5 177 684 A (HARKER PATRICK T ET AL) 
5 January 1993 (1993-01-05) 
figures 2-5,9-18 
column 4, line 50-70 
column 5, line 1-20 

DE 195 40 834 A (SIEMENS AG) 
7 May 1997 (1997-05-07) 



page 2, 
page 4, 
page 5, 
page 6, 



line 1-70 
line 20-65 
line 1-70 
line 1-20 



[X]' 



or documents are listed In 



m 



Patent family members are listed in an 



• Special calegones of 



Hngdate 

.* document which may throw doubts on priority e 
which Is cited to establish the publication date 
dlalion or other special reason (as speciRetO 

3" document referring to an ( 

to the 



'* later document pubKshed after the international fifing date 
or priorSy date and not in oonfilot with the application but 
cited to understand the principle or theoiy underlying the 



relevance; the claimed Inventfon 
._ novel or cannot be considered to 
step when the document is taken alo 



■&* document member of the same p^nt family 



the actual completion of the 

4 April 2001 



11/04/2001 



European Patent Office, P.B. 5818 Patentlaan 2 
NL-22eO HV Riiswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 



Foim PCT/1SA«10 (second sheet) (July 18 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



PCT/US 01/00080 



C.(Contmuatlon) DOCUMENTS 



TO BE RELEVANT 



th incfication.where appropriate, of thi 



Relevant to claim Nc 



US 5 936 517 A (YEH SHOW-WAY) 
10 August 1999 (1999-08-10) 
figures 2,12-19 
column 1, line 35-55 
column 2, line 60-70 
column 3, line 1-60 



Forni PCT/1SA/21D (conlinuJilon of second sheei) (July 1992) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 



Inlerri lal Application No 

PCI/US 01/00080 



Patent document 
cited In search report 


Publication 
date 


Patent family 
memt)er(s) 


Publication 
date 


US 5177684 A 


05-01-1993 


AU 


644664 


B 


16-12-1993 






AU 


8020091 


A 


25-06-1992 






CA 


2046984 


A 


19-06-1992 



DE 19540834 A 


07-05-1997 


AT 


180074 T 


15-05-1999 






AU 


1717297 A 


22-05-1997 






WO 


9716780 A 


09-05-1997 






DE 


59601907 D 


17-06-1999 






EP 


0859988 A 


26-08-1998 






GR 


3030804 T 


30-11-1999 



US 5936517 A 10-08-1999 NONE 



1 POT/1SW210 (patent family annex) (July 19S2) 



